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FOREWORD 


One  aspect  of  the  weapons  effects  research  at  the  Ballistic  Research 
Laboratory  (BRL)  is  the  study  of  penetration  mechanics  for  shaped 
charges  impacting  military  targets.  In  the  Continuum  Mechanics  Center 
(CMC)  of  the  BRL  this  research  involves  theoretical  and  numerical  model¬ 
ing  studies  pertaining  to  the  fundamental  characteristics  and  behavior 
of  materials  subjected  to  dynamic  loads.  The  product  of  this  research 
is  the  development  and  maintenance  of  large-scale  hydrocodes . for  appli¬ 
cation  to  problems  of  military  interest.  The  subject  of  this  report  is 
the  BRL  version  of  the  HELP  code. 
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I.  INTRODUCTION 


One  of  the  objectives  of  the  Continuum  Mechanics  Center  (CMC)  of 
the  Ballistic  Research  Laboratory  is  the  development  and  maintenance 
of  large-scale  continuum  mechanics  codes  at  the  state-of-the-art  for 
application  to  problems  of  military  interest  in  both  basic  research 
and  experimental  studies.  This  objective  is  approached  at  the  CMC  by 
analyzing,  evaluating  and  documenting  existing  calculational  techniques 
and  numerical  models  in  terms  of  validity,  accuracy  and  efficiency,  and 
where  necessary,  refining  and  upgrading  them,  or  suggesting  new  methods 
and  models.  Since  such  considerations  cannot  always  be  the  concern  of 
the  user  looking  for  the  solution  to  a  specific  problem  or  seeking  in¬ 
formation  to  fit  into  some  overall  scheme,  the  CMC  maintains  a  library 
of  working  continuum  mechanics  codes  available  for  use  at  the  BRL,  along 
with  relevant  documentation. 

The  HELP  (Hydrodynamics  ELastic-Plastic)  code1  or  one  of  its  evolu¬ 
tionary  antecedents  has  been  used  for  some  time  as  a  calculational  pre¬ 
dictive  tool  in  research  involving  intermediate  ballistics  (e.g.  muzzle 
blast),  warhead  mechanisms  (e.g.  shaped  charge  jet  formation)  and  ter¬ 
minal  ballistics  (e.g.  impact  and  penetration)  with  varying  degrees  of 
success.  Effective  use  of  the  code  has  been  hampered  at  times  by  the 
memory  size  and  speed  of  the  computers  available  for  running  HELP  cal¬ 
culations,  and  more  seriously,  by  questions  concerning  the  accuracy  of 
some  of  its  predictive  models. 

Furthermore,  many  problems  with  understanding  and  using  a  large 
complex  code  arise  from  a  lack  of  complete  documentation.  When  a  code 
evolves,  such  as  HELP  has,  over  a  period  of  twenty  years,  much  infor¬ 
mation  about  changes  and  corrections  vanishes  and  becomes  part  of  the 
"folklore"  of  the  code.  A  special  effort  will  therefore  be  made  at  the 
CMC  to  document  any  modifications  made  to  such  codes. 

This  report  describes  the  current  effort  at  the  CMC  to  test  and 
evaluate  the  capabilities  and  potential  of  the  HELP  code  as  a  predictive 
tool  in  shaped  charge  studies.  In  Section  II  we  briefly  describe  the 
HELP  model,  the  implementation  of  the  code  on  a  CDC  7600  computer  and 
its  application  to  shaped  charge  design.  For  readers  unfamiliar  with 
shaped  charge  theory  we  recommend  the  work  of  Birkhoff  et  al.  i-n  reference 
2.  Section  III  describes  calculations  made  with  the  code  and  their  com¬ 
parison  with  experiments  and  other  numerical  results .  Section  IV  con¬ 
tains  a  discussion  of  experimental  temperature  measurements  for  shaped 


1  L.  J.  Hageman,  et  al.,  "HELP,  A  Multi-material  Eulerian  Program  for 
Compressible  Fluid  and  Elastic-Plastic  Flows  in  Two-Dimensions  and  Time," 
Systems,  Science  and  Software  Report  SSS-R-75-2654,  July  1975. 

2  G.  Birkhoff,  et  al.,  "Explosion  with  Lined  Cavities,"  J.  Appl.  Phys. 
19,  563,  1948. 
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charges  and  the  difficulties  encountered  with  HELP'S  internal  energy 
calculation.  Conclusions  and  recommendations  for  upgrading  the  code 
are  presented  in  Section  V. 


II.  THE  HELP  CODE 
A.  Brief  Description  of  the  Model 

The  HELP  code  is  a  two-dimensional  Eulerian  code  capable  of 
describing  unsteady  multi-material  interactions  and  treating  material 
strength  as  an  elastic-plastic  phenomenon.  An  important  feature  of  the 
code  is  its  utilization  of  massless  tracer  particles  to  define  material 
interfaces  and  free  surfaces.  These  particles  are  initially  placed 
along  the  surface  boundaries  and  moved  thereafter  with  a  local  material 
velocity.  This  Lagrangian  style  description  of  moving  surfaces  provides 
an  accurate  definition  of  the  material  interfaces  and  free  surfaces 
without  sacrificing  the  ease  of  treating  extreme  distortion  by  an 
Eulerian  model. 

The  HELP  algorithm  numerically  models  the  solution  of  the  con¬ 
servation  equations  governing  the  motion  and  interaction  of  continuous 
media.  In  integral  form  over  a  control  volume  V,  the  equations  are: 


I?  =  'jCPUinidS’ 

lfJ(oujdV  iUjnidS-  3 

ft-4pETdV  ‘fs° ijVidS  -  _^H*iVidS- 


(conservation  of  mass) 


=  1,2,  (conservation  of 
momentum) 


(conservation  of  total 
energy) 


(1) 


(2) 


(3) 


where  the  summation  convention  is  used.  Here  p,  u. ,  a...  E„,  t  denote 

l  lj  T 

the  density,  ith  component  of  the  velocity,  (i,j)  component  of  the 
stress  tensor,  specific  total  energy,  and  time,  respectively.  The  sur¬ 
face  of  the  control  volume  is  denoted  by  S  and  the  ith  component  of  the 

outward  normal  to  it  is  denoted  by  n. .  The  stress  tensor  is  taken  as 

i 


the  sum  of  the  hydrostatic  stress  -6.^P  (6^.  denotes  the  Kronecker  delta 
and  P,  the  pressure)  and  the  stress  deviator  tensor  ^ ;  that  is, 


S.  . 
il 


(4) 


a. . 

il 


6.  .P. 

il 


12 


The  numerical  treatment  of  pure  cells  (containing  only  one 
material)  and  of  mixed  cells  (containing  more  than  one  material)  can 
differ  substantially  in  the  HELP  algorithm.  For  simplicity  and  clarity, 
the  treatment  of  pure  cells  only  will  be  described  below.  The  control 
volume  V  is  taken  to  be  the  computational  cell.  The  time  step  At  is 
determined  using  a  Courant  condition.  Equations  (1)  -  (3)  are  finite 
differenced  over  the  time  step  At  and  can  be  rewritten,  using  equation 


(4),  as: 

SPHASE 

HPHASE 

TPHASE 

A  (m) 

-At 

/  pu.n.dS, 

ys  1  1 

(5) 

A (mu.)  =  At  /  S..n.dS  -  At  / 

1  Js  1  s 

Pn.dS  -At  / 

•  J  Js 

(pUjU^n-dS,  j  = 

1,2, (6) 

A(mET)  =  At  f  S^u-n^S  -  At  J 

1  Pu.n.dS  -  At  / 
c  1  1 

(puiET)nidS. 

(7) 

In  solving  these  equations  the  values  of  the  cell-centered  mass,  velocity 
and  specific  total  energy  at  a  new  time  level  are  found  from  the  values 
at  a  previous  time  level.  This  is  accomplished  in  three  stages  by  deter¬ 
mining  the  time  rate  of  change  of  the  mass,  momentum  and  total  energy 
due  to  the  effects  of:  (1)  the  deviator  stress  (involving  those  terms 
with  the  variable  S. .),  (2)  the  pressure  (involving  those  terms  with  the 

variable  P),  and  (3)  the  transport  of  material  (involving  the  convection 
terms).  These  phases  are  appropriately  named  SPHASE,  HPHASE  AND  TPHASE, 
respectively.  During  each  time  step,  each  value  of  the  mass,  momentum 
and  total  energy  is  updated  sequentially  by  each  phase  in  the  order 
listed  and  each  phase  uses  the  previously  updated  values  as  its  initial 
values.  The  surface  integrals  in  equations  (5)  -  (7)  are  evaluated  using 
the  value  of  the  integrands  at  the  surfaces  of  the  computational  cells. 

Besides  the  dependent  variables  already  listed,  the  specific 
internal  energy  and  the  pressure  need  to  be  determined.  The  former  is 
evaluated  at  any  phase  by  subtracting  from  the  specific  total  energy  the 
specific  kinetic  energy,  that  is, 

Ej  =  Et  -  0.5  (u![  +  ti ^)  • 

Equation  (8)  is  evaluated  at  each  cell  and  at  the  end  of  each  phase. 

The  pressure  is  calculated  prior  to  SPHASE  via  the  algebraic  equation 
of  state  which  has  the  functional  form  P  =  P(p,  E^.) . 

Input  to  the  code  consists  of  an  initial  geometric  configura¬ 
tion,  initial  values  and/or  explosive  loading  parameters,  relevant 
material  properties  and  desired  editing  features.  Output  includes 
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pressure,  density,  material  velocity  components,  specific  internal 
energy  and  deviatoric  stress  components  as  functions  of  space  and  time. 

The  latest  documented  version^  of  the  HELP  code  contains  sev¬ 
eral  features  designed  to  enhance  the  code's  ability  to  model  high 
explosive-metal  interactions  for  complex  geometric  configurations  such 
as  found  in  shaped  charge  designs.  The  geometry  package  is  capable  of 
describing  any  configuration  decomposable  into  straight  line  segments, 
arcs  of  circles,  or  arcs  of  ellipses.  A  high  explosive  detonation 
model  can  account  for  multiple  explosives  and  multiple  detonation  points 
as  Well  as  for  the  presence  of  wave  shaping  devices.  Sliplines  may  be 
employed  at  interfaces.  These  are  especially  useful  at  high  explosive- 
metal  interfaces  where  there  is  a  large  difference  in  material  density. 
The  sliplines  allow  different  materials  in  the  interface  cells  to  devel¬ 
op  velocities  commensurate  with  their  material  properties,  rather  than 
requiring  all  material  in  the  interface  cells  to  move  with  the  same 
velocity. 

B.  Installation  of  HELP  on  the  CDC  7600  Computer 

Before  the  acquisition  of  the  BRL  CYBER  7600,  locally  available 
computers  were  BRLESC  and  UNIVAC  1108.  The  advantages  in  using  the  CDC 
7600  to  perform  calculations  with  a  code  of  HELP'S  complexity  arise 
essentially  from  machine  speed.  A  CDC  7600  possessing  about  160K  octal, 
60-bit  words  of  small  core  and  1.2  million  octal  bits  of  large  core  has 
potential  for  running  problems  with  very  fine  spatial  gridding  over 
long  real  time  periods. 

A  UNIVAC  1108  card-image  tape  version  of  the  code  was  obtained 
from  Systems,  Science  and  Software  to  create  the  current  BRL  7600  ver¬ 
sion  of  the  HELP  code.  The  entire  source  tape  was  initially  transferred 

to  an  UPDATE3  Program  Library,  and  then  transferred  one  subroutine  at 
a  time  via  *ADDFILE  and  *DECK  directives  to  a  new  program  library.  Thus 
each  subroutine  is  a  "deck"  on  the  "program  library",  and  may  be  altered 
and  compiled  separately  from  the  entire  program. 

In  addition  to  the  original  subroutine  decks  (Table  A-l),  three 
new  deck  IDENTS  (Table  A-2),  and  several  "correction"  IDENTS  (Table  A-3) 
were  generated  during  the  installation.  A  discussion  of  each  IDENT 
appears  in  Appendix  A. 

Since  the  BRL's  CYBER  76  computer  system  was  not  initially 
available,  a  telephone  dial-up  connection  had  to  be  instituted  to  a  dis¬ 
tant  host  site.  In  order  to  obtain  the  graphics  display  (both  static 
and  dynamic)  often  necessary  to  interpreting  hydrocode  results,  a  method 
was  devised  to  send  data  telephonically  to  an  in-house  UNITECH  terminal 
driven  tape  drive.  This  data  tape  was  then  used  with  the  BRL  CALCOMP 


3.  Control  Data  Corporation,  UPDATE  Reference  Manual,  7600  Computer 
System,  Arden  Hills,  Minnesota,  1975. 
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plotting  facilities  to  produce  static  results,  or  with  an  interactive 
CRT  to  look  at  dynamic  description  of  phenomena  such  as  liner  collapse 
and  jet  formation  and  break  up.  Figure  1  indicates  schematically  the 
transmission  procedure,  and  the  necessary  programs  for  this  procedure 
are  available  from  CMC  on  request. 

In  the  process  of  checking  out  the  CDC  7600  version  of  HELP, 
two  sample  calculations  were  made  to  compare  running  times  of  the  CDC 
with  those  of  the  BRLESC  and  UNIVAC  1108  computers  for  this  code.  In 
the  first  problem,  the  calculation  of  the  liner  collapse  of  a  81.3  mm 
(called  3.2  inch  precision  shaped  charge)  conical  liner,  the  CDC  7600 
ran  at  a  rate  of  two  seconds  per  time  cycle,  while  the  comparable  BRLESC 
calculation  took  one  minute  per  time  cycle.  Thus,  for  this  problem,  the 
CDC  7600  would  require  about  one  minute  to  perform  a  HELP  calculation 
requiring  a  half  hour  on  BRLESC.  The  second  calculation,  an  iron  rod 
impacting  on  aluminum,  ran  seven  times  faster  on  the  CDC  7600  than  on 
UNIVAC  1108.  It  should  be  noted  that  the  81.3  mm  shaped  charge  calcu¬ 
lation  was  not  run  on  the  UNIVAC  1108  because  this  particular  machine's 
memory  was  too  small  to  permit  a  grid  structure  of  the  size  used  in  the 
calculation  on  the  CDC  7600  and  BRLESC. 

While  running  the  HELP  code  on  the  CDC  7600  several  alterations 
were  made  to  correct  the  coding,  improve  the  code's  input  format  and 
permit  further  output  options.  These  modifications  are  documented  in 
Appendices  A  (see  Table  A-3)  and  B  of  this  report.  Of  particular  note 
is  a  modification  made  to  subroutine  UVCALC.  This  subroutine  is  called 
after  the  HPHASE  calculations  and  updates  the  material  velocities  in 
multi -material  cells  containing  sliplines  by  solving  four  simultaneous 
linear  equations.  (See  reference  1,  pp.  5. 2-5. 3).  When  the  angle  made 
by  the  slipline  relative  to  the  positive  x-axis  is  0°  or  180°,  the 

solution4  of  the  four  equations  is  wrong.  The  IDENT  CORUVC  corrects 
this  error. 

In  addition  to  the  above  mentioned  documentation,  a  set  of  sam¬ 
ple  input  and  corresponding  output  is  provided  in  Appendix  C  to  aid 
users  of  the  BRL  version  of  the  code  (to  be  known  as  BRLHELP) .  Spe¬ 
cifically  the  1/0  corresponds  to  calculations  of  (1)  a  43  mm  conical 
shaped  charge,  (2)  a  hemispherical  shaped  charge,  and  (3)  a  copper  wedge 
impacting  on  a  reflective  wall. 

C.  Application  to  Shaped  Charge  Design 

HELP  has  the  potential  to  simulate  much  of  shaped  charge  beha¬ 
vior,  and  thus  provide  a  needed  numerical  tool  able  to  perform  optimiza¬ 
tion  and  prefabrication  studies  for  materials  and  geometries  proposed 
in  warhead  design.  The  code  has  the  capability  of  modeling  the  collapse 
of  the  liner  and  the  formation  and  development  of  the  jet  and  the  slug. 


4.  A.  R.  Kiwan,  Private  Communication. 
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Figure  1.  J1ELP  Output  Procedure. 


The  code  provides  information  which  can  be  used  to  determine  (1)  jet  tip 
velocity,  (2)  collapse  velocity  and  collapse  angle,  (3)  liner  position 
as  a  function  of  time,  (4)  stretching  rate  of  the  jet,  and  (5)  strength 
effects . 


To  date,  the  code  has  been  used  with  success  in  predicting  jet 

tip  velocities  and  liner  collapse  for  a  variety  of  geometric  designs.5,6 
However,  it  overestimates  the  internal  energy  in  the  collapsed  portion 
of  the  conical  shaped  charge  liner  materials.  In  an  earlier  version  of 

g 

the  code  (BRLSC  -  Ballistic  Research  Laboratory  Shaped  Charge)  ,  this 

overestimation  led  to  unreasonably  low  estimates  of  jet  densities. 

Because  of  this,  an  artificial  bound  on  the  internal  energy  in  the  liner 

was  used  in  the  subroutine  for  the  metal  equation  of  state  (Tillotson 

9  •  . 

Equation  of  State)  in  an  attempt  to  maintain  a  more  appropriate  jet 

density.  Without  this  bound,  the  jet  tip  would  exhibit  internal  ener¬ 
gies  so  high  and  densities  so  low  for  some  liners  as  to  imply  vaporiza¬ 
tion.* *  Until  the  code  is  able  to  predict  internal  energy  accurately, 
internal  energy  dependent  quantities  such  as  temperature,  and  strength 
effects  cannot  be  calculated  with  much  confidence.  A  cause  of  this 
problem  will  be  identified  and  solutions  suggested  in  Section  IV  of  this 
report . 


5. 

ter 

BRL 


J.  T.  Harrison,  "Comparison  Between  the  Eulerian  Hydrodynamic  Compu- 
Code  (BRLSC)  and  Experimental  Collapse  of  a  Shaped  Liner," 

Memorandum  Report  ARBRL-MR- 02841,  June  1978.  ( 


6.  Janet  Lacetera,  "Study  of  Liner  Collapse 
Various  Hemispherical  Shaped-Charge  Systems," 


and  Jet  Formation  for 
BRL  Draft  Report. 


7. 

and 

No. 


A.  Kiwan,  and  A.  Arbuckle,  "Study  of  Liner  Collapse,  Jet  Formation 
Characteristics  from  Impulsive  Shaped-Charge  Systems,  BRL  Report 
2028,  November  1977.  (AD  #A051342) 


o  m  t  Gittines  "BRLSC:  An  Advanced  Eulerian  Code  for  Predicting 
Shaped  Charges.  Volume  I."  BRL  Contract  Report  279.  December  1975. 

(AD  #A023962)  . 

9.  J.  H.  Tillotson,  "Metallic  Equation  of  State  for  Hypervelocity 
Impacts,"  General  Atomic  Report  No.  GS-3216,  July  1962. 


*It  should  be  noted  that  the  artificial  energy  bound  may  be  applied  to 
any  portion  of  a  metal  in  a  shaped  charge  problem,  whether  or  not  that 
portion  of  the  metal  constitutes  part  of  the  liner.  This  is  because 
the  equation  of  state  subroutine  EQST  checks  only  that  the  parameter 
NLINER  (see  Appendix  B)  is  nonzero  before  applying  the  bound  to  the 
appropriate  material;  and  NLINER  is  a  parameter  correlated  with  a  mate¬ 
rial  as  opposed  to  a  material  package. 
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III.  HELP  CALCULATIONS 


Over  the  years  possible  causes  of  the  HELP  code's  shortcomings  have 
been  suggested  as  resulting  from  incorrect  or  incomplete  physical  model¬ 
ing,  incorrect  coding  or  a  combination  of  these  factors.  An  added  com¬ 
plication  heretofore  was  that  HELP  had  been  run  on  computers  with  either 
limited  memory  or  speed.  This  necessitated  the  use  of  various  equiva- 
lencing  and  overlaying  techniques  and  relatively  coarse  grid  structures 
which  might  have  obscured  or  compounded  the  code’s  actual  weaknesses. 
Using  the  CDC  7600  computer  with  its  relatively  high  speed  and  memory 
capacity  diminished  these  complications.  Therefore  to  substantiate 
HELP'S  capabilities  and  determine  its  weaknesses,  several  lengthy  calcu¬ 
lations  made  on  the  CDC  7600  are  reviewed.  The  input  and  selected  output 
from  these  studies  are  presented  in  Appendix  C. 

A.  43  mm  Conical  Shaped  Charge 

The  first  study  using  BRLHELP  modeled  the  liner  collapde  and 
jet  formation  for  a  43  mm  hollow,  conical  copper  liner  driven  by  a  right 
circular  cylinder  of  unconfined  Composition  B  explosive  (Figure  2).  This 

shaped  charge  design  was  chosen  since  both  experimental  and  numerical 
results  existed  for  comparison:  Boyd  Taylor  of  BRL  had  produced  a  series 
of  early-time  flash  radiographs  which  permitted  a  comparison  between 
HELP  tracer  particle  plots  and  photographs  of  the  early  stages  of  the 
liner's  collapse,  as  well  as  estimates  of  the  relative  jet  tip  velocity 
(the  velocity  difference  between  the  tip  of  the  jet  and  the  tail  of  the 
slug).  John  Harrison,  using  the  BRLSC  code,  had  also  performed  an  exten¬ 
sive  comparison  with  this  experiment.^  He  had  found  that  the  code-sim¬ 
ulated  collapse  agreed  quite  well  with  experiment  and  that  the  code  pre¬ 
dicted  that  after  approximately  11.5  ys,  the  jet  tip  reaches  a  constant 
relative  velocity  of  6.54  mm/ys  which  matched  the  experimentally  meas¬ 
ured  value.  Figure  3  illustrates  these  results.  He  also  noted  the 
overly  low  jet  density  and  extremely  high  internal  energy  predicted  by 
BRLSC. 


The  calculations  made  with  the  HELP  code  also  match  the  kine¬ 
matic  results  obtained  experimentally  as  indicated  in  Figure  4.  The 
relative  jet  tip  velocity  calculated  was  6.31  mm/ys  which  differs  from 
Harrison  and  Taylor's  results  by  about  3.3  percent.  This  difference 
results  from  the  bound  on  the  internal  energy  used  in  the  equation  of 
state  which  enforces  a  higher  density  in  the  jet  material  and  conse¬ 
quently  a  lower  velocity.  It  should  be  noted  that  both  calculations 
were  made  without  using  the  strength  options  available  in  each  of  the 
codes  (i.e.  the  subroutine  SPHASE  in  the  HELP  code  and  PH3  in  BRLSC  were 
not  used) . 
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COMPOSITION  B 


Figure  2.  Line  Drawing  of  the  43  mm  Conical  Shaped  Charge. 
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Figure  3.  Comparison  Between  the  BRLSC  Code  Results  and  Experimentally  Measured  Velocity 
Differences  (Jet  Tip  Minus  Slug  Tail).  (Taken  from  Reference  5.) 
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igure  4.  Comparison  of  BRLSC  and  BRLHELP  Code  Results  for  the  43  mm  Conical  Shaped  Charge. 


B.  Hemispherical  Shaped  Charge 

Another  configuration  often  used  in  shaped  charge  design  is 

the  hemisphere.  Since  one-dimensional  codes  such  as  BASC  cannot  pro¬ 
vide  accurate  liner  collapse  and  jet  formation  information  for  designs 
based  on  this  shape,  it  is  very  important  to  have  a  two-dimensional  code 
which  can  serve  as  a  reliable  modeling  tool. 

To  determine  HELP'S  ability  to  model  hemispherical  shaped 
charges,  a  study  was  made  using  a  scaled  version  of  a  five  inch  copper 
hemisphere  for  which  series  of  flash  radiographs  (made  by  Seymour  Kronman 
of  BRL)  existed  to  provide  comparison.  In  this  case  the  copper  liner 
possessed  an  inner  diameter  of  35.814  mm,  a  thickness  of  1.143  mm,  and 
was  driven  by  an  unconfined  right  circular  cylinder  of  Composition  B 
explosive  where  the  head-height  was  19.05  mm,  and  the  charge  was  sub- 
calibered  by  2.54  mm.  The  calculation  was  carried  out  without  using  the 
strength  option  in  HELP.  The  liner  configuration  is  pictured  in  Figure 
5. 


Comparison  of  the  HELP  calculated  tracer  plots  with  photographs 
from  the  flash  radiograph  series  showed  good  agreement  in  size  and  shape 
as  the  jet  developed.  (See  Figures  6  and  7.)  A  systematic  time  differ¬ 
ence  between  the  flash  radiographs  and  the  HELP  results  was  observed. 

The  HELP  calculation  took  three  microseconds  longer  for  the  jet  to  reach 
the  same  state  of  development  as  that  exhibited  in  the  radiographs.  This 
difference  is  due  to  the  fact  that  the  experimental  time  measurements  do 
not  begin  until  a  wire  stretched  across  the  open  edge  of  the  hemisphere 
is  broken  by  the  oncoming  detonation  wave,  while  HELP  starts  counting 

time  from  the  onset  of  detonation  in  the  burn  routine. 

The  relative  jet  tip  velocity  calculated  by  HELP  was  3.06  mm/ys 
while  experimentally  a  value  of  3.4  mm/y  sec  was  measured.  Thus  HELP'S  6 
prediction  was  10  percent  low,  but  well  within  experimental  error  bounds. 

C.  Copper  Wedge 

The  previous  examples  of  HELP  calculations  model  three  areas  of 
shaped  charge  studies,  the  detonation  of  the  explosive,  the  liner  col¬ 
lapse,  and  the  jet  formation.  These  two-dimensional  unsteady  calcula¬ 
tions  are  lengthy  and  complicated.  A  solid  wedge  (conical  or  planar] 
impinging  on  a  reflective  wall  simulates  the  liner  collapse  and  the  jet 

formation2  ,  and  on  occasion,  is  more  suited  for  primary  studies.  Fur¬ 
thermore,  it  does  provide  a  vehicle  to  examine  a  somewhat  different  be¬ 
havior  with  respect  to  the  energy  calculation.  In  problems  described 


TOT  J.  T.  Harrison,  "BASC:  A  Simplified  Analytical  Shaped  Charge  Code," 
BRL  Draft  Report. 
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5.  Line  drawing  of  the  Hemispherical  Shaped  Charge. 
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Figure  6.  Photographs  of  the  Kronman  Flash  Radiographs  Illustrate 
the  Collapse  and  Jet  Formation  for  a  Hemispherical  Liner. 
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t  =  16  ps  t  = 18  ms 

Figure  7.  HELP  Tracer  Particle  Plot  of  the  Liner  Collapse  and  Jet  Formation  for  the 
Hemispherical  Shaped  Charge  Used  in  the  Kronman  Experiment. 


in  Sections  III-A  and  III-B,  the  detonation  of  the  explosive  releases 
internal  energy  generating  a  driving  force  which  accelerates  the  liner 
from  rest  while  internal  energy  is  simultaneously  evolved.  In  the  wedge 
impact  problem  the  wedge  initially  has  all  kinetic  energy  and  the  inter¬ 
nal  energy  is  evolved  from  the  impact  and  deformation  of  the  material. 

At  early  times,  the  physics  of  the  problem  becomes  similar  to  the  shaped 
charge  calculation',  however  the  initial  energy  partition  lends  an  inde¬ 
pendent  check  on  the  source  of  the  internal  energy  problem. 

The  planar  wedge  impact  problem  consists  .of  a  pure  copper  wedge 
at  a  60°  angle  of  obliquity  impacting  a  wall  (the  ordinate)  with  a  velo¬ 
city  perpendicular  to  the  wedge's  length  and  of  magnitude  2.5043  mm/ys 
(see  Figure  8) .  We  note  that  the  direction  of  the  velocity  is  different 
from  that  given  in  reference  2.  However,  a  Galilean  transform  at  the 
intersection  of  the  wedge  and  the  wall  gives  the  appropriate  direction 
to  the  velocity.  The  impact  as  depicted  in  the  calculation  was  per¬ 
formed  without  strength  effects. 


IV.  INTERNAL  ENERGY  PROBLEM 

The  material  properties  of  a  shaped  charge  jet  and  of  a  kinetic 
energy  penetrator  are  important  quantities  to  the  designers  and  engineers 
working  in  the  field  of  penetration  mechanics.  The  accurate  prediction 
of  these  properties  depends  directly  on  the  accurate  calculation  of  the 
internal  energy  of  the  material.  However,  at  the  present  time  the  inter¬ 
nal  energies  predicted  by  the  HELP  computer  code  can  be  unphysically  high 
and  imply  vaporization  in  cases  where  experimental  evidence  is  to  the 
contrary.  In  this  section  we  discuss:  (1)  the  results  of  experimental 
temperature  measurements  for  copper  and  tin- lead  eutectic  conical  shaped 

charge  jets11,  and  (2)  the  solution  to  the  HELP  internal  energy  problem. 

A.  Comparison  with  Experiments 

The  temperature  measurement  of  copper  conical  shaped  charge 
jets  has  been  performed  by  Von  Holle  and  Trimble  at  BRL,  using  two- 

1  2 

color  infrared  radiometry  .  They  made  temperature  measurements  at  the 
jet  tip  and  as  a  function  of  distance  behind  the  jet  tip  which  indicated 
that  the  copper  jets  were  in  the  solid  state. 


11.  W.  G.  Von  Holle  and  J.  J.  Trimble,  "Temperature  Measurement  of 
Copper  and  Eutectic  Metal  Shaped  Charge  Jets,"  BRL  Report  2004,  August 
1977.  (AD  #B021338L) . 

12.  W.  G.  Von  Holle  and  J.  J.  Trimble,  "Residual  Temperature  of  Shocked 
Solids  by  Two-Band  Infrared  Radiometry,"  BRL  Memorandum  Report  2624, 

May  1976.  (AD  #B011450L) 
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Figure  8.  Initial  Conditions  for  Copper  Wedge  Impact  Problem. 
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Measurements  were  made  for  jets  formed  with  both  Composition  B 
and  octol  charges.  Jet  tip  temperatures  for  the  Composition  B  charges 
averaged  714. 15K  with  a  standard  deviation  of  80K,  and  for  the  octol 
charges  averaged  794. 15K  with  a  standard  deviation  of  40  K.  These  tem- 

13 

peratures  are  all  well  below  the  melting  point  (1356. 6K)  of  copper 

Temperatures  behind  the  jet  tip  were  not  as  well  correlated  as 
those  at  the  jet  tip.  In  the  case  of  standard  rounded-apex  liners, 
observations  were  hampered  by  large  amounts  of  debris  near  the  jet. 

Jets  from  specially  machined  conical  liners  were  free  of  debris  behind 
the  jet  tip;  however,  the  scatter  in  the  data  was  quite  large  and,  in 
fact,  most  of  the  temperature  measurements  were  outside  the  calibration 
range  (500- 1100K)  of  the  experiments.  This  fact  is  illustrated  in  Fig¬ 
ure  9  which  was  adapted  from  a  similar  figure  in  reference  11.  However, 
all  measured  jet  temperatures  were  well  below  the  melting  point  of  cop¬ 
per,  and  we  feel  that  comparisons  with  hydrodynamic  theory  and  similar 
experiments  with  tin- lead  eutectic  jets  give  supporting  evidence  that 
the  copper  jets  were  in  the  solid  state. 

The  reliability  of  the  shaped-charge  temperature  measurements 

12 

is  supported  by  the  success  of  Von  Holle's  residual  temperature  mea¬ 
surements  on  copper  plates.  Experimentally  determined  temperatures  of 
copper  plates,  plotted  as  a  function  of  estimated  pressure,  are  compared 

14 

(Figure  10)  with  a  theoretical  curve  due  to  McQueen  and  Marsh  ,  based  on 
bulk  hydrodynamic  theory.  The  theoretical  curve  gives  consistently 
lower  temperatures,  indicating  that  the  experimentally  determined  copper 
plate  temperatures  were  probably  not  underestimates,  and  may  be  consi¬ 
dered  reasonably  accurate.  Furthermore,  we  may  assume  that  those  tem¬ 
perature  measurements  for  the  conical  liners  are  also  reasonably  accu¬ 
rate,  and  assert,  as  did  Von  Holle,  that  the  jets  are  solid  at  least 
near  the  tip. 

Further  evidence  that  the  copper  jets  were  in  the  solid  state 

comes  from  comparisons  with  tin-lead  eutectic  jets'*''*'  which  were  con¬ 
sidered  to  be  in  the  liquid  state.  Framing  camera  and  x-ray  pictures 
of  these  jets  support  a  model  of  the  jets  as  a  condensed  jet  preceded 
by  a  cloud  of  much  less  dense  materials.  A  rather  long  rise  time  (4 
ysec)  of  the  infrared  signal  from  the  jet  lead  to  the  assumption  that 
the  leading  cloud  was  translucent  and  behaved  as  a  volume  radiator. 

Under  the  assumptions  of  a  radiation  model,  consistent  with  the  framing 
camera  and  x-ray  observations,  the  derived  temperatures  indicate  that 


13.  D.  R.  Stull  and  H.  Prophet,  Project  Directors,  JANAF  Thermochemical 
Tables,  2d  ed.,  June  1971. 

14.  R.  G.  McQueen  and  S.  P.  Marsh,  "Equation  of  State  for  Nineteen 
Metallic  Elements  from  Shock-Wave  Measurements  to  Two  Megabars," 

J.  Appl.  Phys.  31^,  1253,  1960. 
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Figure  9.  Temperature  vs  Distance  Behind  Tip  for  Octol  Charges.  (Adapted  from  Reference  11 
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Figure  10.  Plot  of  Copper  Residual  Temperature  Versus  Hugoniot  Pressure 
Exhibiting  Comparison  of  Experimental  with  Theoretical  Data. 
(Adapted  from  Reference  12.) 
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the  cloud  material  and  the  following  condensed  jet  were  both  in  the 
liquid  state.  The  radiation  model  applicable  to  the  copper  jets  was  not 
as  complex,  a  rapid  rise  time  for  the  infrared  signals  indicating  radia¬ 
tion  from  an  opaque  object.  Comparisons  of  this  model  with  that  appli¬ 
cable  to  the  eutectic  jet  observations  indicate  that  the  copper  was  not 
in  the  liquid  state. 

For  the  interested  reader.  Von  Holle  cited  further  evidence 
that  the  copper  jets  were  in  the  solid  state^  and  that  the  eutectic 
jets  were  in  the  liquid  state^. 

In  the  copper  wedge  we  note  that  while  internal  energies  above 
that  required  for  incipient  vaporization  are  evolved,  the  impact  velo¬ 
city  (1.71  km/s)  is  well  below  those  determined  experimentally  which 
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would  induce  this  behavior.  Experiments  by  Weihrauch  ,  with  copper 
projectiles  impacting  copper  targets,  indicate  that  in  the  range  of 
impact  velocities  of  1-3  km/s,  fluid- like  behavior  is  manifested,  with 
surface  temperatures  in  excess  of  melt  (1436-1473K  for  impact  velocities 
on  the  order  of  1.15  km/s).  However,  phase  changes  leading  to  vaporiza¬ 
tion  should  not  occur  until  the  range  of  impact  velocities  of  3-12  km/s 
is  reached. 

Internal  energies  for  both  the  43  mm  conical  shaped  charge  and 
copper  wedge  impact  problems  generated  by  the  HELP  code  are  high  enough 
to  imply  that  portions  of  the  collapsed  liner  and  deformed  wedge  are  in 
the  vapor  state  (see  Figures  11  and  12).  On  the  basis  of  the  experimen¬ 
tal  evidence  such  energies  are  considered  unphysical.  The  artificial 
bound  on  the  internal  energy  used  in  the  subroutine  for  the  metal  equa¬ 
tion  of  state,  for  the  shaped-charge  problem,  does  maintain  greater  jet 
densities  than  would  be  obtained  without  the  bound.  However,  internal 
energy  dependent  variables  such  as  temperature  and  material  strength 
remain  unphysical.  In  the  case  of  the  hemispherical  shaped  charge,  the 
calculated  internal  energies  do  not  approach  incipient  vaporization. 

This  does  not  imply  that  the  calculated  internal  energy  is  correct,  but 
only  that  we  have  no  obvious  indication  of  error  (see  Figure  13) . 


15.  R.  Karpp  and  J.  Simon,  "An  Estimate  of  the  Strength  of  a  Copper 
Shared  Charge  Jet  and  the  Effect  of  Strength  on  the  Breakup  of  a 
Stretching  Jet,"  BRL  Report  1893,  June  1976.  (AD  #B012141L) 

16.  G.  Hauver,  "Residual  Temperature  Measurements  on  Shock  Compressed 
Metals,"  Bull.  Am.  Phys.  Soc.  20_,  19,  1975. 

17.  G.  Weihrauch,  "The  Behavior  of  Copper  Pins  Impacting  Materials  with 
Velocities  Between  50  m/s  and  1650  m/s,"  Dissertation  toward  the  academic 
degree  of  Doctor-Engineer  in  Mechanical  Engineering,  University  of 
Karlsruhe,  February  1971. 
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Figure  11.  Specific  Internal  Energy  of  Copper  for  43  mm  Conical  Shaped  Charge  Computed  Along  the 
Axis  of  Symmetry  at  15ys  Using  BRLHELP. 


4.0 


u 

u 

UJ 

a, 

to 


1  11.0mm 


SLUG  END 


JET  TIP 


Figure  12.  Specific  Internal  Energy  for  Copper  Wedge  Impact  Computed 
Along  the  Axis  at  15ys  Using  BRLHELP 
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Figure  13.  Specific  Internal  Energy  for  a  Hemispherical  Shaped  Charge 
Computed  Along  the  Axis  of  Symmetry  at  16ys  Using  BRLHELP. 
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B.  Numerical  Analysis  of  HELP!s  Energy  Formulation 

It  is  the  purpose  of  this  section  to  identify  the  sources  of 
error  in  the  internal  energy  calculation  and  examine  possible  solutions. 
The  algorithm  by  which  the  specific  internal  energy  is  calculated  in 
pure  cells  is  outlined  in  Section  II-A.  We  recall  that  the  specific 
internal  energy  at  a  given  time  level  is  computed  as  the  difference 
between  the  specific  total  energy  and  the  specific  kinetic  energy.  The 
former  is  updated  using  a  finite  difference  approximation  to  the  conser¬ 
vation  of  total  energy  equation  (3)  and  the  latter  is  obtained  by  an 
algebraic  manipulation  of  the  updated  values  of  the  mass  and  momentum 
derived  from  finite  difference  approximations  to  the  conservation  of 
mass  and  momentum  equations  (1)  and  (2),  respectively.  All  of  the  above 
quantities  are  calculated  at  the  center  of  the  computational  cell. 

The  numerical  treatment  of  mixed  cells  is  far  more  complicated 
and  in  some  instances  substantially  different  from  the  pure  cell  treat¬ 
ment.  As  an  example,  the  cell  pressure  in  a  multimaterial  cell  is 
obtained  by  an  iterative  procedure  whereby  for  a  constant  specific  inter¬ 
nal  energy  value  of  each  material  the  densities  are  varied  until  the 
pressures  converge  to  the  same  value.  Furthermore,  in  the  TPHASE  cal¬ 
culation,  the  properties  of  each  material  in  a  slipline  cell  are  updated 
using  previously  determined  individual  material  properties  within  the 
cell  rather  than  the  cell-centered  properties.  These  specific  treatments, 
although  important  and  indispensable  to  the  correct  running  of  the  code, 
make  an  accurate  and  complete  analysis  of  the  numerical  treatment  of 
mixed  cells  unwieldy.  Consequently,  the  following  discussion  will  address 
only  the  pure  cell  algorithm. 

Two  possible  sources  of  error  in  the  internal  energy  calcula¬ 
tion  are  (1)  the  functional  evaluation  via  equation  (8)  of  the  internal 
energy,  and  (2)  the  subtraction  of  two  large  quantities  to  obtain  a  small 
quantity.  Since  the  former  error  tends  to  be  small  in  actual  calcula¬ 
tions,  and  since  the  latter  error  would  be  restricted  to  at  most  a  small 
portion  of  a  calculation  (for  the  class  of  problems  considered),  neither 
error  would  seem  to  be  a  major  contributor  to  a  gross  error  in  the  inter¬ 
nal  energy. 

Another  possible  source  of  error  is  the  numerical  algorithm  used 
to  compute  the  total  energy  or  the  kinetic  energy  or  both.  After  a  de¬ 
tailed  study  of  the  total  energy  algorithm,  no  error  or  discrepancy  was 
found  in  the  manner  in  which  it  was  updated  at  each  phase.  However,  this 
was  not  the  case  with  the  kinetic  energy  calculation.  To  demonstrate  the 
problem  associated  with  the  kinetic  energy  algorithm,  we  will  initially 
limit  our  discussion  to  that  portion  of  the  HELP  calculation  that  contri¬ 
butes  most  significantly  to  the  artificial  value  of  the  internal  energy- 
TPHASE.  The  following  statements  will  be  proven:  (1)  the  approximations 
used  in  TPHASE  to  update  the  mass  and  momentum  of  a  cell  are  first  order 
accurate  approximations*  to  the  TPHASE  portion  of  the  partial  differential 

*A  finite  difference  approximation  to  a  differential  equation  is  said  to 
be  first  order  accurate  if  the  formal  truncation  error  is  first  order. 
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equations  governing  conservation  of  mass  and  momentum;  (2)  the  calcula¬ 
tion  of  the  corresponding  kinetic  energy  includes  terms  which  are  of  the 
order  of  the  truncation  error,  (3)  the  inclusion  of  these  first  order 
terms  produces  smaller  kinetic  energies  than  that  of  a  strictly  first 
order  approximation  and,  thus,  larger  internal  energies  via  equation  (8) • 

Before  establishing  the  order  of  the  mass  and imomentum  approxi¬ 
mations,  we  will  outline  the  mass  and  momentum  approximations  in  TPHASE. 
In  order  to  simplify  the  discussion  we  consider  the  one-dimensional 
Cartesian  formulation.  Recall  that  TPHASE  considers  only  the  transport 
of  material  from  cell  to  cell  and  does  not  consider  any  terms  associated 
with  the  pressure  or  deviator  stress.  Thus,  in  modeling  this  phase  the 
conservation  of  mass  and  momentum  in  one  dimension  is 


|f  +  |-  (pu)  =  0,  (9) 

,  1;.  [(pu)u]  .  0.  (10) 

Consider  three  adjacent  cells  i-1,  i  and  i+1  of  equal  length  Ax.  The 
center  of  the  ith  cell  is  at  x  =  (i-^)Ax.  Let  the  velocity  be  positive; 
that  is,  the  material  is  moving  to  the  right.  See  Figure  14.  The  inter¬ 
face  between  cell  i-1  and  i  is  denoted  by  i  and  that  between  i  and  i+1 
by  r.  Integrating  equations  (9)  and  (10)  over  cell  i,  we  obtain 


9_ 

9t 


/  pdx 
*/cell  i 


(pu)r 


(pu)£  =  0, 


[(pu)u]^  =  °> 


(ID 

(12) 


where  the  subscripts  r  and  £  denote  the  evaluation  of  the  bracketed 
quantities  at  the  interfaces  r  =  iAx  and  l  =  (i-l)Ax.  The  integrals  in 
equations  (11)  and  (12)  represent  the  mass  nr  =  p^Ax  and  the  momentum 

miui  =  piuiAx  cel1  respectively.  Equations  (11)  and  (12)  are  the 

one-dimensional  TPHASE  portion  of  equations  (1)  and  (2),  respectively. 
The  HELP  code  approximates  equations  (11)  and  (12)  by  equations 


m1?4^  -  nu  +  6m  -  6m  =  0, 
i  i  r  S, 


n+1  n+1  n  n  ,  n. 

m.  u.  -  m.u.  +  u. 6m 
li  ii  i  r 


Ui-16\ 


=  0, 


(13) 

(14) 
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Figure  14.  Geometry  Associated  with  the  TPHASE  Calculation  in  One 
Dimension  with  a  Positive  Velocity. 


where  superscript  n  denotes  post-HPHASE  time  level  and  n+1  the  post- 
TPHASE  time  level.  During  a  timestep  At,  an  amount  of  mass  6m  and  momen¬ 
tum  u6m  associated  with  cells  i-1  and  i  will  be  transported  to  cells  i 
and  i+1,  respectively,  since  the  velocity  is  positive.  Furthermore, 

the  HELP  algorithm  approximates  the  convected  masses  6m  ,  6m  by 

i  r 


6,\  =  pi-iVt> 

(15) 

6m  =  p?u  At, 
r  i  r 

(16) 

where  u^  and  ur  are  the  corresponding  transport  velocities  and  a  unit 

cross  sectional  area  is  assumed.  We  note  that  the  specific  momentum  of 
each  6m  is  based  on  its  cell-centered  velocity. 

We  shall  now  prove  that  the  HELP  approximations  (13)  -  (16) 
form  a  first  order  accurate  approximation  to  equations  (9)  -  (10)  within 
TPHASE.  We  take  for  example  the  conservation  of  mass.  Using  the  rela¬ 
tionship  m  =  pAx  and  equations  (15)  and  (16),  we  rewrite  (13)  as 
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p  .u  - 
l  r 
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P  •  ,u. 

'l-l  Z 


At 


Ax 


=  0. 


(17) 


For  this  simple  analysis  we  shall  use  the  average  velocity  between  cell 
i-1  and  i  for  u^  and  i  and  i+1  for  u^.  Although  the  HELP  code  uses  a 

more  sophisticated  approximation  to  u  and  u  ,  the  qualitative  results 

are  the  same.  We  expand  p?+1,  p"^,  u^+1  in  a  Taylor  series  with 

a  remainder  about  the  ith  cell  center  and  the  nth  time  level  and  obtain 
to  second  order 
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where  0,  a,  6,  y  are  between  zero  and  one.  Substituting  the  above 
expansions  into  equation  (17),  we  have 
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From  equation  (19),  we  see  that  the  finite  difference  approximation  (17) 
differs  from  the  partial  differential  equation  (9)  by  terms  multiplied 
by  Ax  or  At.  Therefore,  equation  (17)  is  a  first  order  accurate  approxi¬ 
mation  to  equation  (9).  The  formal  truncation  error  is  O(Ax)  and  O(At). 
However,  the  approximation  (17)  is  equivalent  to  the  HELP  approximation 
(13)  and,  consequently,  the  HELP  approximation  (13)  is  said  to  be  a  first 
order  approximation  to  the  conservation  of  mass  equation  (9) .  When  equa¬ 
tion  (13)  is  used  to  approximate  equation  (9),  we  observe  from  equation 
(19)  that  all  O(Ax),  O(At)  and  higher  order  terms  are  excluded.  We  shall 
see  that  this  is  not  true  in  the  kinetic  energy  calculation.  Similar 
comments  and  observations  to  those  on  the  conservation  of  mass  approxi¬ 
mation  hold  true  for  the  conservation  of  momentum  approximation. 

The  governing  equation  for  the  kinetic  energy  (obtained  from  an 
algebraic  manipulation  of  equations  (9)  and  (10))  is 

—Cpe)  +  (pu)e  =  0,  (20) 


where  e  =  u  /2.  To  calculate  the  specific  kinetic  energy,  the  HELP  algo¬ 
rithm  takes  the  already  computed  momentum  nu divides  it  by  the 

n+i 

value  of  itk  ,  squares  it  and  multiplies  the  result  by  0.5.  To  determine 

the  accuracy  of  the  HELP  approximation  to  equation  (20) ,  we  calculate  the 
specific  kinetic  energy  in  terms  of  the  HELP  approximations  (13)  and  (14) 
From  equation  (14)  we  have 


_  c  /  n+l\  „  c/  n+l\  f  n  n  n  n  1 

0.5  u.  )  =  0.5  m.  )  m.u.  -  u.6m  +  u.  ,6m 

\  i  /  V1/  L11  ir  1-1  ju 

We  expand  the  bracketed  term,  use  equation  (13)  and  obtain 


(21) 
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The  value  obtained  from  equation  (22)  for  the  specific  kinetic  energy 
is  identical  to  the  value  the  HELP  algorithm  calculates.  We  now  perform 
a  similar  analysis  on  equation  (22)  to  that  performed  on  equation  (13). 

2 

By  using  equations  (15),  (16),  m  =  pAx  and  e  =  u  / 2,  we  rewrite  equation 
(22)  in  a  more  convenient  form: 
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We  now  expand  the  quantities  p?+1,  pi  l*  ui  l*  ui+l'  +  1  and  (p©)^  j 

in  Taylor  series  about  the  cell  center  i  and  time  level  n,  substitute 
them  into  equation  (23)  and  obtain 
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The  constants  a,  3,  y,  and  a  are  between  zero  and  one  and  are  associated 

with  Taylor  expansions  of  (pe)V_^,  u?_^,  ui+l’  anc*  *-pe^i  »  respectively. 

In  order  to  retrieve  equation  (20)  from  the  right-hand  side  of  equation 
(24),  0(Ax)  and  0(At)  terms  must  be  neglected.  Thus,  equation  (23)  is  a 
first  order  approximation  to  equation  (20) .  We  note  that  the  terms 

0. 5|pu(3u/8x)  Ax  and  -0.  sj^pu^  (8u/8x) At  in  equation  (24)  are  the 

lowest  terms  derived  from 


n  n 

Vi-Ipi 

9  n+1 

2P  i 


(Ax  -  u^At) 


(25) 


and  are  of  the  order  of  the  truncation  error.  Thus,  expression  (25) 
which  is  included  in  HELP'S  kinetic  energy  approximation  is  of  the 
order  of  the  truncation  error  (0(Ax),  0(At)).  If  one  would  neglect 
expression  (25)  in  equation  (23),  the  results  would  still  be  a  first 
order  approximation  to  equation  (20) . 


Although  in  the  theoretical  limit  as  the  mesh  approaches  zero 
the  truncation  error  term  (25)  approaches  zero,  in  practice  the  inclusion 
of  it  in  equation  (23)  alters  the  computed  value  of  the  kinetic  energy  of 
each  cell  at  each  cycle.  Some  important  properties  of  expression  (25) 
are:  (1)  The  factor  (Ax  -  u  At)  is  positive  and  the  sign  of  the  whole 

term  depends  on  the  transport  velocity  u^.  Consequently,  it  can  be 

easily  shown  that  it  decreases  the  value  of  the  kinetic  energy  at  the 
new  time  level.  This  decrease  is  converted  to  an  increase  in  the  inter¬ 
nal  energy  via  equation  (8).  (2)  In  regions  of  steep  velocity  gradients, 
the  magnitude  of  expression  (25)  increases  significantly  because  of  the 
squared  factor.  In  a  conical  shaped  charge  calculation,  this  term  can 
decrease  a  cell's  kinetic  energy  by  as  much  as  15  percent  in  one  time- 
step.  During  a  calculation  which  typically  consists  of  hundreds  of 
cycles,  the  effect  on  the  internal  energy  by  this  truncation  error 
expression  in  the  kinetic  energy  calculation  at  each  time  level  is 
cumulative  and  may  be  substantial.  (3)  Expression  (25)  does  not  model 
any  term  in  the  partial  differential  equation  (20)  governing  the  kinetic 
energy.  It  is  solely  a  result  of  the  finite  difference  approximation 
(13)  and  (14).  If  we  were  to  finite  difference  equation  (20)  directly, 
it  would  be 


n+1  n+1  n  n  »  n 
mi  ei  =  miei  + 


T  11 

.  .  -  6m  e. . 
l-l  r  i 


(26) 


Comparing  equations  (22)  and  (26),  we  find  that  the  last  term  of  equation 
(22),  which  is  equivalent  to  expression  (25),  is  excluded. 
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Furthermore,  a'similar  analysis  to  the  one  above  for  the  TPHASE 
portion  of  the  calculation  can  also  be  done  for  the  SPHASE  and  HPHASE 
portions  and  a  similar  result  obtained.  The  terms  corresponding  to 
expression  (25)  in  HPHASE  and  SPHASE  depend  on  the  density  and  the  gra¬ 
dient  squared  of  the  pressure  and  normal  stress,  respectively.  However, 
they  are  only  O(At)  and  would  increase  the  kinetic  energy  (the  internal 
energy  would  decrease) .  However,  in  practice  the  value  of  At  is  much 

8  ~2 

less  than  the  value  of'Ax[0(10  )  vs  0(10  )]  .  Consequently,  the  effect 

of  these  extraneous  terms  may  not  be  as  significant  as  in  TPHASE.  For 
example,  in  a  similar  situation  which  gave  15  percent  decrease  to  the 
kinetic  energy  in  TPHASE,  the  corresponding  terms  in  HPHASE  gave  only  a 
0.25  percent  increase. 

We  have  shown  that  the  terms  of  the  order  of  the  truncation 
error  which  are  included  in  the  kinetic  energy  calculation  increase  the 
computed  values  of  the  internal  energy.  This  algorithm  mechanism  is  the 
cause  of  the  unphysical ly  high  internal  energy  values  seen  in  conical 
shaped  charge  simulations.  A  detailed  analysis  of  the  two-dimensional 
algorithm  and  the  effects  of  the  truncation  error  terms  on  the  accuracy 
of  the  internal  energy  calculation  is  given  in  reference  18.  We  note 
that  these  truncation  error  terms  are  direct  consequences  of  the  method 
in  which  the  kinetic  energy  is  updated  via  the  finite  difference  equa¬ 
tions  for  mass  and  momentum  in  HELP.  If  one  retains  the  total  energy 
formulation  in  HELP  and  neglects  these  terms  in  the  calculation  of  the 
internal  energy,  the  present  kinetic  energy  algorithm  must  be  altered. 
The  other  alternative  is  calculating  the  internal  energy  directly  from 
the  partial  differential  equation  governing  internal  energy. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  preceding  sections  have  sought  to  briefly  describe  the  BRL  CDC 
version  of  the  HELP  code,  and  to  identify  some  of  its  strong  and  weak 
points  in  relation  to  shaped-charge  behavior.  We  showed  in  Sections 
IIItA  and  III-B  that  velocities  andshapes  of  developing  jets  are  accu¬ 
rately  computed.  However,  we  know  from  Section  IV-A  that  the  internal 
energy  computed  by  HELP  is  unphysically  high  and  from  Section  IV-B  the 
cause.  With  these  facts  in  mind,  we  outline  a  program  to  upgrade  the 
code's  performance  while  retaining  its  best  features. 

The  cause  of  the  unphysical  internal  energies  has  been  identified 
as  the  inclusion  of  terms  in  the  computation  of  the  kinetic  energy, 
which  are  of  the  order  of  the  truncation  error,  in  the  HELP  algorithm. 


18.  J.  A.  Schmitt,  "A  New  Internal  Energy  Calculation  for  the  HELP  Code, 
and  Its  Implications  to  Conical  Shaped  Charge  Simulations,"  BRL  Report 
ARBRL-TR-02168,  June  1979,  and  to  appear  in  J.  Comp.  Phys.'  (AD  #A072785) 
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Terms  of  this  order  are  only  included  directly  in  the  calculation  of  the 
kinetic  energy  and  thus,  implicity,  via  equation  (8),  in  the  calculation 
of  the  internal  energy.  We  recommend  that  another  method  which  neglects 

1 8 

these  terms  be  incorporated  into  the  HELP  code  .  Two  possible  approaches 
are:  (1)  a  direct  internal  energy  formulation  based  on  the  partial  dif¬ 

ferential  equation  for  internal  energy,  or  (2)  retention  of  the  current 
total  energy  formulation  but  deriving  computed  values  of  the  kinetic  en¬ 
ergy  from  the  partial  differential  equation  for  kinetic  energy,  and  not 
from  the  updated  values  of  mass  and  momentum.  In  either  case  we  expect  a 
more  accurately  calculated  value  of  the  internal  energy,  and  hence,  a 
more  useful  HELP  code. 

Accurate  calculation  of  the  internal  energy  in  metallic  liners  will 
make  viable  the  employment  of  a  more  sophisticated  metal  equation  of 
state  in  the  HELP  code.  Such  an  equation  of  state  is  the  improved 
19  20 

(BRL-GRAY)  version  of  the  GRAY  equation  of  the  state.  It  is  a  multi¬ 
phase  (solid- liquid-vapor)  equation  of  state  which  treats  mixed  phases 
in  explicit  detail,  and  is  therefore  particularly  suited  to  the  treat¬ 
ment  of  jet  formation  and  breakup. 

The  BRL-GRAY  EOS  will  be  installed  in  the  HELP  code  and  comparison 
of  this  EOS  with  the  Tillotson  EOS  will  be  made  for  thermodynamic  regimes 
not  well  described  by  the  Tillotson  EOS.  Thus,  the  requirements  for  an 
EOS  such  as  BRL-GRAY  for  such  calculations  may  be  established. 

We  will  also  be  making  a  modification  to  the  HELP  code  that  will 
eliminate  most  of  the  pressure  iteration  problems  that  we  have  encoun¬ 
tered  in  our  running  of  the  code,  and  that  have  plagued  other  users  of 
the  code. 

Ultimately  we  hope  to  provide  a  numerical  capability  with  HELP  that 
can  be  used  reliably  in  prefabrication  studies  for  proposed  liner  mate¬ 
rials  and  gometries  and  generally  provide  theoretical  support  for  experi¬ 
mental  efforts. 


19.  Joseph  Lacetera,  "BRLGRAY:  The  Ballistic  Research  Laboratory  Version 
of  the  GRAY  Equation  of  State  with  Mixed-Phase  Data  for  Al,  Fe,  Cu,  Ta 
and  U,”  BRL  Draft  Report. 

20.  E.  M.  Royce,  ”GRAY,  A  Three-Phase  Equation  of  State  for  Metals,” 
Lawrence  Livermore  Laboratory  Report  No.  UCRL-51121,  September  1971. 
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MAIN 

EDIT 

FRACS 

NEIVRHO 

PTSAV 

TPHASE 

CARDS 

I  DENT 
HELPCOM 

HELPEQ 

HELP75 

I  DENT 
BRL76 

CHINPUT 

CHKTIM 

CML77 

COMMOD 

CORUVC 

ECARDS 


APPENDIX  A.  UPDATE  Modifications 


Table  A- I:  Subroutine  DECK  Names 


ADDTCR 

CALFRC 

CDT 

COMPRSN 

COMPRS 

DETIME 

DMADJ 

ENCHCK 

ENDMV 

EOUT 

EQST 

ERROR 

FILGRD 

FLGSET 

HPHASE 

INFACE 

INPUT 

MAP 

MOVTCR 

NEWFLG 

NEWMIX 

P LG ADD 

PLGALF 

PLGGEN 

PLGMAS 

PLGTCR 

PLGVOL 

PLUGUV 

REZONE 

RNDOFF 

SETUP 

SETUPA 

SPHASE 

STRNG 

THETAS 

TSETUP 

UVCALC 

UVMOD 

VDCLOS 

XPTOC 

YPTOC 

ADDENG 

UC 

VC 

VOLFND 

WN 

WT 

XCTOP 

YCTOP 

Table  A-II:  Additional  DECK  Names 
EXPLANATION 


COMMON  deck  containing  blank  COMMON,  labeled  COMMON, 
and  related  DIMENSION  statements. 

COMMON  deck  containing  EQUIVALENCE  statements  and  an 
additional  DIMENSION  statement. 


DECK  name  used  when  originally  loading  source  program 
from  card  images. 


Table  A-II I:  Correction  IDENTs 
EXPLANATION 


Simplifies  a-variable  I/O;  adds  print  suppression  options 
deletes  subroutine  CARDS.  (See  Appendix  B) 

Corrects  FORMAT  error  in  original  coding. 

Enables  code  to  halt  calculation  before  running  out  of 
CPU  time. 

Suppresses  the  repetitious  CDC-FORTRAN,  informative  error 
message  115  by  using  the  CDC  utility  subroutine  SYSTEMC. 
The  diagnostic  is  caused  by  the  argument  of  the  intrinsic 
function  EXP  being  smaller  than  -675.84. 

IDENT  which  was  used  when  installing  HELPCOM  and  HELPEQ 
decks  into  program  library. 

Corrects  algebraic  error  in  the  subroutine  UVCALC.  (See 
Section  II-B) 

Added  EQUIVALENCE  statement  to  subroutine  CARDS  for 
operation  on  CDC  7600. 
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ident  explanation 


ENENCHK 

Provides  printing,  in  each  cycle,  of  kinetic  and  internal 
energy  summs  for  all  the  cells  in  the  grid  for  both 
material  and  cell -centered  values. 

JOHN 

Corrects  original  coding  in  JWL  EOS  in  subroutine  EQST. 

LCHANGE 

3 

Corrections,  made  by  Laura  Hagerman  of  S  and  transmitted 
to  TBD,  dealing  primarily  with  program  errors  and  correc¬ 
tions  to  the  automatic  void  closing  routine. 

MEMEXPND 

Modifies  HELPCOM  to  increase  the  number  of  cells  in  a 
calculation. 

MOD76,  MOD76A 

Used  in  loading  subroutines  from  HELP75  deck  to  new 
program  library. 

NOEDPRT 

Provides  options,  via  NAMELIST  input,  for  suppressing 
certain  printout  from  subroutine  EDIT . 

PRT 

Added  print  statements  to  HPHASE  and  CMPRSN  to  overcome 
error  in  CDC-FORTRAN  4.3  compiler.  This  error  does  not 
exist  in  the  CDC-FORTRAN  4.6  compiler. 

PRTPLT 

Copies  variables  to  a  deck  file  for  subsequent  plotting. 

PRTPLTL 

Changes  PRT  prints  to  WRITE  on  FILE  9. 

RID 

Removes  improperly  labeled  sections  of  coding. 
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APPENDIX  B.  Modified  Input 


Instructions  for  generating  HELP  problems,  using  the  UNIVAC  1108 
version  of  the  code,  are  given  in  Section  7.2  of  reference  1,  Here  we 
indicate  input  differences  between  the  BRL7600  and  UNIVAC  1108  versions 
of  the  code. 

In  the  HELP  code  input  parameters  located  in  the  first  150  words 
of  blank  common  are  called  Z-block  variables.  They  are  so-called  because 
they  follow  Z(l),  the  first  word  of  blank  common,  and  therefore  can  be 
referenced  as  Z(I),  where  I  can  take  any  value  between  1  and  150.  In 
the  UNIVAC  1108  version  these  parameters  are  read  by  a  subroutine  called 
CARDS.  In  the  BRL7600  version  this  subroutine  has  been  eliminated  and 
the  Z-block  variables  are  read  using  the  NAMELIST  option. 

A  set  of  edit  control  parameters  are  first  read  using  the  NAMELIST 
OPTNS.  The  Z- variables  are  then  read  using  the  NAMELISTS  START  and  RUN. 
These  are  described  below*. 

NAMELIST  OPTNS:  (Edit  control  parameters.  It  should  be  noted 
that  the  first  three  of  these  parameters  can  be  used  to  suppress  diag¬ 
nostic  print  and  should  be  used  with  caution.) 

IFLGST:  =  0,  default; 

=  1,  suppresses  message  "SUM  OF  FRACS  NOT  EQUAL  TO 
TOTAL  AREA  I,  J,  TAU,  XDY2PI  =  ". 

IPMADJ:  =  0,  default; 

=  1,  suppresses  message  "EVALUATION  OF  MATERIAL  N 
I,  J,  N,  M,  ML,  MB,  TFLUX,  XMASS  (N,M) ,  SAMMP  (N,M) , 
SAMPY  (N,M)" . 

ITPHASE:  =  0,  default; 

=  1,  suppresses  message  "MASS  EVAPORATED  DUE  TO 
ROUND-OFF  IN  TPHASE  I,  J,  MFLAG,  AMX,  T,  E" . 

NPRTOP:  =  0,  default; 

=  1,  suppresses  edit  on  sliplines; 

=  2,  suppresses  long  edit  on  last  cycle; 

=  3,  combines  effects  of  (1)  and  (2). 

TIMMAX:  The  maximum  time  in  CPU  seconds  (decimal)  that  the 

code  will  use  in  the  calculations. 

For  NPRTOP  =  2  or  3  TIMMAX  should  be  about  4  seconds 
less  than  the  JOB  CARD  time;  for  NPRTOP  =  1  or  0 
TIMMAX  should  be  about  10  seconds  less  than  the 
JOB  CARD  time. 

NAMELIST  START 


PK(1):  Problem  number. 
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Restart  cycle  number. 


PK(2) : 

PK(3):  =  0,  for  initial  run; 

=  -1,  for  restart  without  short  edit  print  of  restart 
cycle; 

=  -2,  for  restart  with  short  edit  print  of  restart 
cycle.  (see  p.  7.21  of  reference  1) 

NAMELIST  RUN  (Contains  all  Z -block  variables  not  contained 
in  NAMELIST  START). 


Input  Z- Block  Variables 


Variable  Default 

Name  Value 


BBAR  . 5 


CRATIO 


CVIS 


CYCMX  2 


CYCPH3  1 


Definition 


A  constant  used  in  the  calculation  of  the 
local  sound  speed  of  materials  other  than 
ideal  gases  and  high  explosives.  C  =  Cq  + 

BBAR  •  \[\P\  . 

The  compression  of  the  material  in  adjacent 
cells  is  computed.  If  the  ratio  of  these 
cells'  compressions  is  greater  than  CRATIO, 
their  pressures  are  inverse  compression 
weighted  and  their  velocities  are  compres¬ 
sion  weighted  to  define  cell  boundary 
values  in  HPHASE. 

A  flag  that  describes  the  condition  of  the 
bottom  grid  boundary.  If  CVIS  =  0,  the 
bottom  boundary  will  be  reflective.  If 
CVIS  =  -1,  the  bottom  boundary  will  be 
transmittive. 

The  number  of  passes  through  subroutine 
INFACE  each  cycle  to  minimize  transport 
noise  near  interfaces. 

The  number  of  passes  through  subroutine 
SPHASE  (strength  phase)  each  cycle.  When 
the  calculation  is  purely  hydrodynamic  (no 
strength  effects),  set  CYCPH3  =  -1. 


SO 


Variable 

Name 

DMIN 


DTMIN 

EMIN 

EMOB 

FINAL 

GAMMA 

IEXTX 

ICSTOP 


Default 

Value 


10 


-11 


Definition 


The  maximum  relative  error  in  the  energy 
sum  that  the  user  wants  to  tolerate.  The 
error  is  computed  as  follows: 

(KMAX  \ 

T  E,  -  ETH  )  /ETH 
k=2  k  / 

where  E^.  is  total  energy  of  cell  k  and 

ETH  is  theoretical  energy  in  the  grid  -- 
computed  in  SETUP  and  updated  in  HPHASE, 
TPHASE,  SPHASE,  ADDENG,  REZONE,  and  RNDOFF . 

The  minimum  value  for  the  timestep.  The 
program  gives  an  error  exit  if  CDT  calcu¬ 
lates  a  At  <  DTMIN. 

The  minimum  value  of  specific  internal 
energy  to  be  used  in  the  ideal  gas  equa¬ 
tion  of  state. 

This  variable  was  read  by  CARDS  as  a  dummy 
end  card  of  a  setup  deck  after  all  other 
cards  were  read.  This  card  is  not  used  in 
BRLHELP  input. 


.4  The  final  value  of  the  stability  fraction 

used  in  determining  the  timestep.  See 
STAB . 

y  in  P  =  (y-1)  E  for  material  #20,  an 
ideal  gas. 

A  rezone  flag.  The  grid  will  be  rezoned 
in  the  x-direction  only  when  IEXTX  =  1. 
Omit  this  card  if  the  grid  is  not  to  be 
rezoned,  or  if  it  is  to  be  rezoned  in  the 
y-direction  only. 

Gives  cycle  at  which  execution  will  stop 
if  the  user  wants  to  specify  a  cycle, 
rather  than  a  value  of  T  (time),  on  which 
to  stop.  If  stopping  on  time,  omit  this 
card. 
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Variable 

Name 


Default 

Value 


Definition 


IGM 

When  IGM  =  0  code  uses  cylindrical  coordi¬ 
nates.  When  IGM  =  1,  code  uses  plane  coor¬ 
dinates. 

IMAX 

The  number  of  columns  in  the  grid.  If 
planning  to  rezone  the  grid  in  the  x-direc- 
tion,  IMAX  must  be  an  even  number.  IMAX 
also  must  be  at  least  3. 

INTER 

A  special  editing  flag  for  debugging: 

=  2,  for  debug  prints  in  CDT 
=  3,  for  debug  prints  by  EDIT  after 
HPHASE  and  SPHASE 
=  5,  for  debug  prints  in  TPHASE 
=  7,  for  debug  prints  in  SPHASE 
=  11,  for  debug  prints  in  UVMOD 
=  13,  for  debug  prints  in  UVCALC. 

IPCYCL 

The  number  of  cycles  between  EDIT  prints. 

If  printing  on  time  intervals,  omit  this 
card. 

IPLGBT 

The  bottom-most  row  of  the  plugging  region 
of  the  target.  (Usually  the  first  row  of 
the  target.)  Omit  this  card  when  not  using 
the  plugging  option. 

IPLGRT 

The  right-most  column  of  the  plugging 
region  of  the  target.  (2-3  columns  to  the 
right  of  the  expected  location  of  the 
plug’s  vertical  edge.)  Omit  this  card 
when  not  using  the  plugging  option. 

IPLGTP 

The  top-most  row  of  the  plugging  region 
of  the  target.  (Usually  the  top  row  of 
the  target.)  Omit  this  card  when  not 
using  the  plugging  option. 

IPR 

35  Maximum  number  of  iterations  CDT  will  per¬ 

form  when  attempting  to  equilibrate  the 
pressures  of  materials  in  a  multi-material 
cell.  If  the  pressures  are  not  within  an 
epsilon  (PRCNT)  of  the  average  pressure 
after  iterations,  an  error  exit  occurs. 
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Variable 

Name 


Default 

Value 


Definition 


11 

The  number  of  columns  in  the  grid  that  have 
nonzero  velocities  or  energy  plus  2.  11 

and  12  define  the  "active"  grid.  (Many 
quantities  are  computed  only  for  cells 
inside  the  active  grid.) 

12 

The  number  of  rows  in  the  grid  that  have 
nonzero  velocities  or  energy  plus  2. 

JEXTY 

A  rezone  flag.  The  grid  will  be  rezoned 
in  the  y-direction  only  when  JEXTY  =  1. 

Omit  this  card  if  the  grid  is  not  to  be 
rezoned,  or  if  it  is  to  be  rezoned  in  the 
x-direction  only. 

JMAX 

The  number  of  rows  in  the  grid.  If  plan¬ 
ning  to  rezone  in  the  y-direction,  JMAX 
must  be  an  even  number.  JMAX  also  must  be 
at  least  3. 

KUNITR 

7  The  name  of  the  file  INPUT  will  read  from. 

KUNITW 

7  The  name  of  the  file  EDIT  and  SETUP  will 

write  on. 

LVISC 

A  linear  artificial  viscosity  term  is  added 
to  cell  boundary  pressures  if  LVISC  =  1. 

MAPS 

When  MAPS  >  0,  part  of  the  EDIT  print  is 
a  set  of  symbolic  maps  of  the  compression 
(or  density),  pressure,  radial  and  axial 
velocities,  and  specific  internal  energy 
of  the  cells  in  the  active  grid. 

MAPS  =  1,  gives  a  compression  map. 

MAPS  =  2,  gives  a  density  map. 

MINX 

MAXX 

MINY 

MAXY 

These  parameters  specify  the  left  and  right 
columns  and  the  bottom  and  top  rows,  respec 
tively,  of  the  region  in  which  ADDTCR  will 
add  tracer  particles.  These  variables  must 
be  defined  when  NADD  /  0. 

NADD 

A  flag  for  automatically  adding  material 
tracer  particles  in  the  region  specified 
by  MINX,  MAXX,  MINY,  MAXY.  If  NADD  =  10 
ADDTCR  is  called  every  cycle  which  is  a 
multiple  of  10,  (e.g.,  cycle  20,  30,  ... 
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Variable 

Name 


Default 

Value 


Definition 


250,  260,  etc.).  NOTE:  NTRACR  (=Z(72)) 
must  also  be  defined  before  material 
tracers  can  be  added  by  ADDTCR. 

NDUMP7 

Gives  frequency  of  restart  tape  dumps 
relative  to  EDIT  prints.  (Program  dumps 
only  when  EDIT  prints.)  e.g. ,  if  NDUMP7  = 
5,  a  restart  dump  will  be  made  every  5th 
time  EDIT  prints. 

NFRELP 

Gives  frequency  of  "long"  EDIT  prints, 
which  give  velocities,  energy,  compression 
(or  density)  and  stresses  for  all  cells  in 
active  grid.  A  "short"  EDIT  gives  this 
information  only  for  the  first  column  of 
cells,  e.g.,  when  NFRELP  =  2,  every  second 
EDIT  print  is  ’’long,”  when  NFRELP  =  5, 
every  fifth  EDIT  print  is  "long”,  etc. 

NMXCLS 

The  maximum  number  of  interface  cells  to 
exist  in  the  grid  on  any  one  cycle  during 
the  calculation.  This  maximum  should  cor¬ 
respond  to  the  dimensions  of  the  material 
arrays  (XMASS,  SIE,  etc.). 

NLINER 

The  package  number  of  the  shaped  charge 
liner  material  that  forms  the  jet.  Used 
only  when  calculating  the  collapse  of  a 
liner. 

NMAT 

The  number  of  material  packages,  excluding 
the  void  package. 

NODUMP 

When  NODUMP  =  1,  EDIT  will  not  write  any 
restart  dumps.  (This  flag  overrides  the 
NDUMP7  option,  but  does  not  prevent  SETUP 
from  writing  the  cycle  0  dump.) 

NOSLIP 

If  no  sliplines  are  to  be  generated,  set 
NOSLIP  =  1,  which  causes  the  instructions 
and  routines  that  affect  only  slipline 
cells  to  be  skipped. 

NSLD 

The  maximum  number  of  cells  the  user 
expects  will  contain  the  si ip line  on  any 
one  cycle.  This  maximum  should  also  de¬ 
pend  on  the  dimensions  of  the  slipline 
arrays  (MSLD,  etc.). 
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Variable 

Name 


Default 

Value 


Definition 


NTCC 

When  NTCC  >  0,  NTCC  passive  tracer  parti¬ 
cles  are  initially  placed  in  the  center 
of  every  fourth  nonempty  cell  and  subse¬ 
quently  moved  through  the  grid  with  the 
material.  These  are  used  only  for  produc¬ 
ing  particle  plots  and  do  not  affect  the 
material  behavior.  When  this  card  is 
omitted  (NTCC  =  0),  cell- centered  tracers 
are  not  computed  unless  the  plugging  option 
is  being  used,  in  which  case  they  are  auto¬ 
matically  generated  in  the  plugging  region 
of  the  target. 

NTPMX 

The  maximum  number  of  tracer  particles  per 
material  package.  This  maximum  should  cor¬ 
respond  to  the  dimensions  of  the  TX  and  TY 
arrays . 

NTRACR 

The  number  of  material  tracers  per  cell 
diagonal.  Used  by  ADDTCR  when  adding 
tracers.  ADDTCR  adds  tracers  only  when 
the  distance  between  two  consecutive  tra¬ 
cers  in  a  specified  region  is  greater  than 
a  cell  diagonal  divided  by  NTRACR.  MINX, 
MAXX,  MINY,  MAXY  specify  the  region  of  the 
grid  considered  by  ADDTCR.  NADD  specifies 
the  frequency  with  which  ADDTCR  is  called. 

NUMREZ 

The  total  number  of  times  the  grid  will  be 
automatically  rezoned.  If  not  automati¬ 
cally  rezoning  the  grid,  omit  this  card. 

NUMSCA 

The  number  of  times  the  cycle  or  time  in¬ 
terval  between  EDIT  prints  is  increased. 
(See  PRLIM) 

NVRTEX 

The  second  index  of  the  void  tracer  parti¬ 
cle  that  is  at  the  vertex  of  the  void  clos¬ 
ing  region,  i.e.,  (TX(NV0ID,  NVRTEX), 

TY (NVOID,  NVRTEX))  are  the  coordinates  of 
the  vertex  point.  Used  only  when  the 
code’s  automatic  void  closing  routine  is 
to  be  activated. 

PK(i) 

Problem  number.  Any  number  from  .0001  to 
99.9999. 
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Variable  Default 

Name  Value 


Definition 


PLGOPT 


PLWMIN 


PMIN  5  x  106 

PRCNT  .001 

PRDELT 

PRFACT 

PRLIM 


If  the  user  is  generating  a  plugging  cal¬ 
culation,  PLGOPT  must  be  set  equal  to  1. 

The  special  plugging  mechanisms  are  acti¬ 
vated  only  if  PLGOPT  =  1.  (See  Chapter  VI, 
Ref.  1,  before  using  the  plugging  option.) 

The  required  minimum  specific  plastic  work 
(ergs/gm)  which  the  material  near  the  top 
of  the  vertical  edge  of  the  plug  must  have 
before  the  plug  edge  is  extended  to  the 
next  row.  Omit  this  card  when  not  using 
the  plugging  option. 

A  pressure  cutoff.  If  |P(k)|  <  PMIN,  then 
then  P(k)  =  0. 

Convergence  criterion  for  the  pressure 
equilibration  of  material  in  a  multi¬ 
material  cell.  All  material  pressures, 

P^,  must  satisfy  the  following: 

I (P  -  P.)/P  I <PRCNT. 

1  av  l  av 1 

The  time  (sec)  between  EDIT  prints.  If 
printing  on  cycle  intervals,  omit  this 
card. 

The  factor  by  which  the  print  interval 
(PRDELT  or  IPCYCL)  and  the  print  limit 
(PRLIM)  are  increased.  (See  PRLIM) 

The  time  or  cycle  at  which  the  EDIT  print 

frequency  is  changed.  Example:  you  wish 

-8  -  7 

to  print  every  10  sec  until  T  =  10 

—  7  —  6 

sec;  and  every  10  sec  until  T  =  10 

-6 

sec,  and  every  10  sec  thereafter.  Set: 

PRDELT  =  10-8 
IPCYCL  =  0. 

PRLIM  =  10"7 
PRFACT  =10. 

NUMSCA  =  2. 

NOTE:  When  PRDELT  is  increased  by  a  fac¬ 
tor  PRLIM  is  also  increased  by  the  same 
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Variable 

Name 


PROB 

REZ 

ROEPS 

SIEMIN 

STAB 

TSTOP 


Default 

Value  _ Definition 


factor  for  the  next  rescaling.  If  you 
want  a  constant  print  interval,  omit 
NUMSCA,  PRLIM,  PRFACT . 

Problem  number.  Is  set  equal  to  PK(1) 
by  the  code. 

The  flag  which  initiates  the  rezoning  of 
the  grid.  If  NUMREZ  >  0,  this  flag  will 
be  set  automatically  by  the  code.  However, 
the  user  can  force  a  rezone  on  any  restart 
cycle  by  setting  REZ  =  1  in  the  restart 
deck.  (IEXTX  and/or  JEXTY  must  also  be 
defined  in  order  for  the  grid  to  be 
rezoned. ) 

10”5  A  round-off  epsilon  used  to  define  cutoffs 

in  the  calculation. 

105  A  cutoff  on  the  total  specific  internal 

energy  increment  of  a  cell  in  TPHASE.  This 
cutoff  prevents  small  numerical  signals 
from  enlarging  the  active  grid. 

10  The  initial  value  of  the  stability  fraction 

used  in  determining  the  time  step.  If 
FINAL  =  0,  STAB  is  constant.  Otherwise, 
its  value  is  doubled  each  cycle  until  it 
reaches  the  value  of  FINAL. 

The  value  of  T  (time)  at  which  the  calcu¬ 
lation  will  stop.  When  stopping  on  a 
specified  cycle  (ICSTOP  >  0) ,  set  TSTOP  = 

0.  TSTOP  must  always  be  included  in  both 
the  setup  deck  and  the  restart  deck. 
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APPENDIX  C.  Sample  Calculations 


This  appendix  contains  listings  of  input  decks  and  corresponding 
output  for  three  sample  problems: 

1.  43mm  conical  shaped  charge; 

2.  Hemispherical  shaped  charge; 

3.  (a)  Copper  wedge  startup;  and, 

(b)  Copper  wedge  restart. 

The  input  deck  listings  are  contained  in  Section  C^I  along  with  an  expla¬ 
nation  of  the  job  control  language  (JCL)  statements  necessary  for  running 
these  problems  on  a  CDC  7600.  The  remaining  sections  of  Appendix  C  deal 
with  the  code  output  for  each  problem. 


SD 


Appendix  C-I:  Input  Deck  Listings 

The  JCL  stream  typically  contains  seven  cards  if  input-output  (10) 
is  performed  on  disk.  These  are  as  follows: 


Card 

1 

"Job  card",  specifies  user  name,  mainframe  station, 
time  limit. 

Card 

2 

Request  for  permanent  file  space  allocation  (for 
code  restart  file,  TAPE  10). 

Card 

3 

GETPF,  which  puts  a  copy  of  the  HELP  program  library 
onto  a  local  file  (A). 

Card 

4 

Full  UPDATE  with  OLDPL  =  A,  NEWPL  =  B. 

Card 

5 

FTN  card. 

Card 

6 

LGO  card. 

Card 

7 

Catalogs  "TAPE  10"  on  permanent  file. 

The  next  cards  shown  in  the  startup  (initial  run)  decks  are:  a 
problem  identifier  card  (Hollerith);  the  NAMELIST  input  described  in 
Appendix  B;  grid  specifications;  material  package  information  including 
initial  density,  velocity,  and  specific  internal  energy  for  each  material 
strength  constants;  tracer  particle  specifications;  slipline  information 
if  used;  and  specifications  of  high  explosive  detonation  points. 

To  run  the  copper  wedge  problem  a  modification  is  necessary  to  the 
basic  geometry.  This  modification  is  shown  in  the  UPDATE  input  stream 
with  the  IDENT  WEDGEOM.  It  shortens  the  YAXIS  array  and  lengthens  the 
XAXIS  array  to  allow  as  large  a  copper  bar  as  possible  within  the  con¬ 
straints  of  the  problem  and  the  total  number  of  cells  allowed  (12001). 

It  should  also  be  noted  that  in  this  problem  no  sliplines  (NOSLIP  =  1) 
are  employed. 

In  the  copper  wedge  restart  one  JCL  card  is  added  to  put  the  restart 
file  onto  a  local  file.  The  code  expects  this  file  to  be  called  TAPE7. 

It  should  also  be  noted  that  the  general  specifications  following  the 
NAMELIST  input  on  initial  runs  is  not  needed  here  as  this  information 
is  already  available  to  the  code  from  the  restart  file.  The  variable 
PK(2)  =  11  specifies  the  restart  cycle  number  and  PK(3)  =  -1  indicates 
that  the  problem  is  a  restart.  ICSTOP  =  14  specifies  that  the  calcula¬ 
tion  should  halt  after  14  cycles. 
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MBR48 , STMFZ , T2000 .  SAMPLE  PROBLEM  2-CU  HEMISPHERE  LACETERA, JANET 
REQUEST (TAPE7,  *PF) 

GETPF (A,BRLHELP, ID=CMC HYDRO) 

UPDATE (F,P^A,N=B) 

FTN ( I=COMPILE , L=0 , PL=55000 ,  I£M=I ) 

RFL (160000) 

LGO. 

CATAIOG (TAPE 7, HEMI , IDCMCHYDRO) 


EXIT. 

SAMPLE  PROBLEM  2  -  COPPER  HEMISPHERE  LOADED  WITH  COMP  B 


$OPTNS  NPRTOP=2 , 

TIMMAX=1000, 

IFDGST=1 , 

IPMADJ=1, 

$START  PK (1)=1, 

PK(2)=0, 

PK(3)=-0$ 

$RUN  TSTOP=l . 5E- 
KUNITW=7 , 
NOSLIP=0 , 

04, 

KUNITR=7 , 
LVISC=1 , 

PRDELT=  2.0E- 

-06, 

NMAT=2, 

NMXCLS=800, 

NTRACR=5 , 

OTPMX=900 

CYCFH3=-1, 

NLINER=1 , 

NSLD=300, 

MAPS=1, 

IPR=100, 

IMAX=60, 

JMAX=200, 

DMIN=100, 

NFRELP=100, 
CVIS=-1 , 

60 

999 

200 

NDUMP7=1, 

IGM=0$ 

.05 

.05 

IPCYCL=0 , 

1 1=3, 

999 


HE  UNCONFINED 
ITPHSE=1$ 


REZ=0. , 


12=53, 


0 

3 

8.94 

21 

1.717 

2.35E+09 

6.95E+10 

0.0 

0.0 

4  1 

300 

3.14159 

1.57080 

0.0 

6.35 

1  1 

50 

1.905 

•  6.35 

-4  1 

300 

1.57080 

3.14159 

0.0 

6.35 

1  2 

100 

0.0 

2.54 

2  2 

200 

2.159 

2.54 

1  2 

50 

2.159 

6.35 

-4  2 

300 

1.57080 

3.14159 

0.0 

6.35 

4  3 

300 

3.14159 

1.57080 

0.0  0.0 

0.0  0.0 

5.50E+10  5.30E+09 

0.0  0.0 

1.905 

1.7907  6.35 

1.7907 

2.159  2.54 

2.159  6.35 

1.905  6.35 

1.905 


0.0 

0.0 


4.55E+11 

0.0 


.9785 

1.0 
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0.0 

6.35 

1.7907 

1 

3 

50 

1.7907 

6.35 

1.905 

6.35 

1 

3 

50 

1.905 

6.35 

2.159 

6.35 

2 

3 

200 

2.159 

6.35 

2.159 

2.54 

-1 

3 

100 

2.159 

100 

2.54 

0.0 

2.54 

1 

0 

1 

0 

300 

0 

0 

2 

0 

301 

0  600 

1 

2 

0.0 

2.54 

2 

2 

1.82 

5.62 

0 

0.0 

2 

.54 

2.159 

6.35 

0.0 

5 

.59 

2.159 

6.35 

*EOI 
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MBR48 , SIMFZ , T77 .  BRLHELP  TEST  WEDGE  1.0  STARTUP 

REQUEST (TAPE 10 , *PF) 

GETPF  (A,  BRLHELP,  IDCMCHYDRO) 

UPDATE (F,P=A,N=B) 

FTN  (I  COMPILE ,  L=0 ,  PL=55000 ,  LCM=I ) 

LGO. 

CATALOG  (TAPE  10 , 60WEGECU,  IDCMC HYDRO) 

*EOR 

*IDENT  08NOV77 
*1  EDIT. 240 

305  CONTINUE 
*D  NOEDPRT.8 
*IDENT  WEDGEOM 
*D  MEMEXPND. 3 

COMMON /MISC/  PQ (5,8) ,XAXIS(5,85) ,YAXIS(5,140) ,VL(4) 

*D  MEMEXPND. 4 

COMMON/XDXYDY/  X00,X(85)  ,Y00,Y(140)  ,DX00,DX(85)  ,DY00,DY(140) 
*D  MEMEXPND. 12 

2  MFLAG (12001) ,TAU(85) ,UL(400) ,PL(400) ,TX(5,900) , 

*D  MEMEXPND. 33 

COMMON/MAP/  PROP (85) ,PR(85) 

*EOR 

COPPER  IMPACT  STUDY  60  DEGREE  WEDGE 
$OPTNS  NPRTOP=3 ,  TIMMAX=  60.,  IFLGST=1,  IPMADJ=1,  ITPHSE=1$ 
$START  PK (1 ) =29 . $ 

$RUN  KUNITW=10, 

NOSLIP=l ,  LVISC=1 , 

NMAT=1 . ,  NMXCLS=800. ,  NTRACR=5 . ,  NTPMX=900,  REZ=0. , 

CYCFH3=-1. , 

MAPS=1. ,  11=85,  12=140,  IMAX=85,  JMAX=140, 

TSTOP=5 . 0E-06 ,  NFRELP=100,  NDUMP7=1,  IPCYCL=50 , 


CVIS= 

-1. 

f 

IGM=1$ 

60 

6 

5 

14 

.1 

.25 

5 

1.0 

999 

10 

10 

100 

10 

.5 

» 

in 

CN 

1 

.25 

10 

.5 

999 

0 

2 

8, 

.94 

0.0 

-2.170E+05 

1. 

,  250E+05 

2.350E+09  6.950E+10  5.500E+10  5.300E+09  4.550E+11  .9785 


3 

0.0 

1  200 
12.0 

19.00 

22.972 

2 

19.0 

1  50 

22.972 

19.0 

24.129 

-3 

19.0 

1  200 
24.129 

0. 

13.159 

3 

0.0 

2  200 
13.159 

19.0 

24.129 

2 

19.0 

2  50 

24.129 

19.0 

22.972 

-3 

19.0 

2  200 
22.972 

0. 

12.0 

100 
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MBR4 8 , STMFZ ,  T77 .  BRLHELP  TEST  WEDGE  2.0  RESTART 

REQUEST (TAPE 10 , *PF) 

GETPF (TAPE7, 60WEGECU,CY=1 , ID=CMC HYDRO) 

GETPF (A, BRLHELP, ID=CMC HYDRO) 

UPDATE (F,P=A,N=B) 

FTN  (I  COMPILE ,  L=0 ,  PL=55000 ,  LCM=I ) 

LGO. 

CATALOG  (TAPE  10 , 60WEGECU ,  IDCMCHYDRO) 

*EOR 

*IDENT  08NOV77 
*1  EDIT. 240 

305  CONTINUE 
*D  NOEDPRT, 8 
*IDENT  WEDGEOM 
*D  MEMEXPND. 3 

COMMQN/MISC/  PQ(5,8) ,XAXIS(5,85) ,YAXIS (5,140) ,VL(4) 

*D  MEMEXPND.  4 

COMMON /XDXYD Y /  X00,X(85) ,Y00,Y(140) ,DX00,DX(85) ,DY00 ,DY (140) 

*D  MEMEXPND. 12 

2  MFLAG (12001) ,TAU(85) ,UL(400) ,PL(400) ,TX(5,900) , 

*D  MEMEXPND. 33 

COMMON /MAP/  PROP (85) ,PR(85) 

*EOR 

COPPER  IMPACT  STUDY  60  DEGREE  WEDGE 
$OPTNS  NPRTOP=3 ,  TIMMAX=  60.,  IFIGST=1,  IPMADJ=1,  ITPHSE=1$ 

$START  PK (1)=29. $ 

$START  PK (1 ) =29 . ,  PK(2)=11.,  PK(3)=-1.$ 

$RUN  KUNITW=10,  KUNITR=7, 

NOSLIP=l ,  LVISC=1 , 

NMAT=1. ,  NMXCLS=800. ,  NTRACR=5 , ,  NTPMX=900,  REZ=0, , 

CYCFH3=-1., 

MAPS=1. ,  11=85,  12=140,  IMAX=85,  JMAX=140, 

TSTOP=0 , 0E-06 ,  NFRELP=100,  NDUMP7=1,  IPCYCL=50,  ICSTOP=14, 

CVIS=-1 , ,  IGM=1$ 

*EOI 
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MBR48 , STMFZ ,  T77 .  43-MM  SHAPED  CHARGE 

REQUEST (TAPE10 , *PF) 

GETPF ( A , BRLHELP , ID=CMC HYDRO) 

UPDATE ( F , P=A , N=B , L=A1234 ) 

FTN ( I=COMPILE , L=0 , PL=55000 , LCM=I ) 

LGO . 

CATALOG (TAPE10 , 43MMCALC , ID=CMCHYDRO) 

43-MM  SHAPED  CHARGE  CALCULATION 

$OPTNS  NPRTOP=2 ,  TIMMAX=60.,  IFLGST=0,  IPMADJ=1,  ITPHSE=0$ 

$START  PK (1 ) =43  .  $ 

$RUN  NOSLIP=0 ,  LVISC=0 ,  DMIN=1.0E03, 

NMAT=2 . ,  NMXCLS=800 . ,  NTRACR=5 . ,  NTPMX=900,  REZ=0., 


MAPS=0 

,  11=  60 

,  12=187,  IMAX=60 , 

JMAX=1 87 , 

T  r  - 

ICSTOP 

=849, 

PRDELT=5.0E-07, 

NLINER=1 , 

NSLD=300 

TSTOP= 

0.0, 

NFRELP=1000 ,  NDUMP7=1, 

IPCYCL=0 

CYC PH  3 

-+1.  , 

KUNITW=10 , 

CVIS=- 

1.,  IGM= 

0$ 

46 

111 

.052  .0572 

.0629 

.0692 

1 

111 

.0761  .0837 

.0921 

.1013 

1 

111 

.1114  .1225 

,1348 

.1483 

1 

1  1 

.15  .1992 

.1993 

999 

187 

999 

0 

.052 

2 

8.9 

0.0  0.0 

0.0 

21 

1.717 

0.0  0.0 

0.0 

2 . 350E+09  6. 950E+10  5.500E+10  5.300E+09  4.550E+11  .9785 

0.0  0.0  0.  0.  0.  1.0 


1 

3.142 

0.000 

100 

1.950 

3.696 

0.254 

1 

0.236 

300 

3.600 

2.296 

8.573 

1 

2.296 

50 

8.573 

2.159 

8.573 

1 

2.159 

300 

8.573 

0.118 

3.649 

1 

1.950 

0.000 

100 

3.142 

3.696 

0.127 

2 

0.000 

100 

1.270 

2.296 

1.270 

2 

2.296 

200 

1.270 

2,296 

8.573 

2 

2.296 

300 

8.573 

0.236 

3.600 

2 

1.950 

0.000 

100 

3.142 

3.696 

0.254 
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4  3 

3.142 
0.000 

100 

1. 

3. 

950 

696 

0.127 

3  3 

0.118 

300 

3. 

649 

2.159 

8.573 

1  3 

2.159 

50 

8. 

573 

2.296 

8.573 

2  3 

2.296 

200 

8. 

573 

2.296 

1.270 

-1  3 

2.296 

100 

1. 

270 

0.000 

1.270 

100 

1  0 

1 

0 

400  0 

0  2 

0 

301 

0  700 

1 

2 

0.0 

1.27 

0 

1 

.  27 

2 

.29646  8.57 
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Appendix  C-II:  43  mm  Shaped  Charge  Output 

The  complete  code  output  is  shown  for  the  43  mm  shaped  charge  cal¬ 
culation,  along  with  one  page  of  UPDATE  output  and  three  pages  of  system 
output.  The  first  two  pages  of  UPDATE  output  were  calls  to  COMDECKS  and 
are  not  shown.  The  third  UPDATE  page  shows  the  correction  IDENT  and 
DECK  lists.  The  system  output  consists  of  the  SCOPE  2  load  map. 

The  first  page  of  code  output  is  the  OPTNS  namelist  followed  by: 
material  package  geometry  and  slidepoint  information;  a  short  edit  on 
the  explosive  detonation*;  Z-variable  values;  initial  conditions  for  each 
material;  material  properties  needed  for  strength  calculations  (note 
that  specific  internal  energy  as  well  as  pressure  are  given  relative  to 
reference  values  at  ambient  conditions);  mesh  specifications;  a  message 

-  8 

that  cell  (1,  25)  has  been  detonated  and  that  t  -  3.2964778  x  10  sec; 

a  cycle  zero  summary;  eight  pages  of  cell  coordinates  for  tracer  parti¬ 
cles;  and  a  summary  of  the  total  energy  in  the  grid  in  terms  of  both 
cell -centered  and  material  values  for  SPHASE  and  HPHASE. 

These  are  followed  by  part  of  a  map  of  the  mixed  and  pure  cells. 

If  the  Z-variable  MAPS  had  not  been  set  to  zero  these  would  also  have 
been  a  complete  set  of  maps  for  the  other  output  variables  (compression, 
pressure,  radial  velocity,  axial  velocity,  and  specific  internal  energy). 

The  remaining  output  is  typical  for  terminal  output  with  NPRTOP  =  2 
and  MAPS  =  0.  Following  the  tracer  particle  edit  is  an  edit  of  output 
parameters  along  the  y-axis  of  the  problem  with  the  y-axis  of  the  prob¬ 
lem  with  the  y-axis  index  varying  from  185  to  1. 

It  should  be  noted  that  while  IFLGST  and  ITPHSE  are  zero  in  this 
problem  the  corresponding  output  does  not  appear  here  because  the  con¬ 
ditions  for  these  outputs  have  not  yet  occurred. 

To  obtain  the  full  print  out,  the  following  coding  must  be  inserted 
in  DETIME  in  place  of  statement  436  CONTINUE: 

436  WRITE (6, 530)  J, (DETIM(K+I) ,  I=1,IMAX). 


*The  full  detonation  edit,  including  detonation  times  for  each  detona¬ 
tion  point,  is  automatically  suppressed  in  this  version  of  BRLHELP  but 
will  appear  as  optional  output  in  subsequent  versions. 
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Appendix  C-III:  Hemispherical  Shaped  Charge  Output 


This  section  includes  only  selected  output;  however,  a  large  number 
of  maps  are  shown  to  indicate  the  detail  one  can  expect  from  this  type 
of  output.  Inspection  of  the  compression  map  for  each  cycle  informs  the 
reader  of  the  location  of  each  material,  as  they  are  separated  by  mixed 
cells  denoted  by  asterisks  (*) . 

For  this  reason  it  was  not  considered  essential  to  reproduce  the 
pure  and  mixed  cell  maps,  which  are  a  part  of  the  normal  output.  All 
other  maps  give  the  output  value  in  the  mixed  cell  as  well  as  the  pure 
cells . 
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TYPE  OF  INITIATION  POInT  1  EXPLOSIVE  PACKAGE  2  INITIATION  POINTS  0.  2.54000E+00  UtLAY  TIME 
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-H-H-M-H-II-H-II-M-H-H-H 

-I-I-i-l-l-l-l-l-l-l-l 
-l-l-l-l-l-I-I-l-l-l-l 
-U-U-U-d-a-l-l-l-1-1-1 
-a-a-d-d-4-a-a-J-u' 

,-A-R-R- J-4-4-4-4-4- 4-4' 


C  0  U  0  6  I  M 
F  E  (  11  f  «  H  H  |  II  J  f  * 
66FC0EFHHHH1  1  IF 
(  (  F  C  t  C  F  (  H  H  H  |l  H  |  H  I 

H  II  «  t  (  E  [  (  F  6  6  6  II  6  H  II  Ml 

h»M66f  FEttF*H  6  6  6IIHH 
I  11  II  H  M  6  6  f  t  t  [  F  F  6  f  F  6  6  B  b  1  6 

J  J  J  1  1  1  II  M  H  6  f  t  E  t  E  F  F  f  f  F  f  F  6  6  6  G 

J  J  J  U  l  I  M  H  H  G  f  f  E  I  E  E  I  f  F  F  f  F  f  f  F  f  G 

djl  1  1  IH66FFEEECEELEECFFFFFF 

JlllMHUfFEEEOOOOEtLEtEEFFFF 
lHM66FFFECEDOnOODUEEEEEEEEF 
6666FFECL0000D000UUUE  EE  E  EEE 
6FFFELCOODCCCECCOOUDOOEEEEE 
FfttEllllOCCCCCCCCCOOllOODDllEE 
EEOODOCCCCBBOCCCCCCOODDUODD 
0UDCCCCVBBB0DD8CCCECD0DDDDD 
CCCC0UUUBU8OBBi1CCCCCCCOUDOD 
CbBbBBAAAAAaBBBBCCCCCCCCCCC 
BbBAAAAAAAAABBBBCEECECCLECC 
AAAAA---AAAACBBBBBCCCCCCCCC 
.  — a-A-A-A-A-A  -AAABbHnBWBUECCCCCC 
A-A-A-A-A-A-A-A  -AAAABBBBBBbBBEECLb 
U-K-b-M-A-A-A-A  -AAAAABBBBBbbbbHbHb 
O-U-B-ri-A-A-A-A  -AAAAAAOnBBiBBBbBDB 
a-b-b-U-A-A- A-A  --AAAaAAABUBB«NHB88 
b-b-b -b- A- A • A- A- A  --AAAAAAAABibBBbBb 
■U-b-U-b- X • A- A - A “A - A- A  ..aAAAAAAAAAUBAA 
b-B-b-b-S-A-A-A-A-A-A-A-A-A  --AAAAAAAAAAA 
•  8-«-B-H-B-U-n-|l-A-A-A-A-A-A-A-A  -  -AAAAAAAAA 
.O-b-B-b-8-U-a-b-b-b-A-A-A-A-A-A-A-A-A  -AAAAAAA 
>C  -C  -  R-H-H-ti-H-H-H-b-b-b  -  8-  A- A  -  A  -  A  -  A- A  .  A  -  A  -  -  -  -  A  A 
•C-C-C -C-C-C-b-H-b-b-U-b-H-A-b-H-A-A- A-A  »A-A - A • A  * A  - A  - A 
.C-C-C-C-C-C-C-C-C-C-B-8-8-U-B-B-B-b*b- A-A-A-A-A-A-A-A. 

.C-C-C-C-C-C-C-C-C-C-C-C-C-l-B-b-a-b-H-b-B-b-H-A-A-A-A 

•a-U-D-U-O-O-C-C-C -C-C-C-C -C-C'C-C -b-b-U-b-B-fc-b-B-b-A 
•o-u-D-u-n-u-u-D-u-n-c-c-c-c-c-c-c-L-c -c-a-H-b-b-n-B-f 
•U-U-O-U-D-U-U-O-U-U-D-D-O-O-C-C-C-C-E -L-C -C -C -C-V-b-B 
•E-E-E-t-E-E-E-O-O-U-U-D-O-n-O-t-O-D-U-C-C-C-C-C-C-C-C 
.E-C-E-t-t-L-E -t-C-E-C-E-C-0-D.U-O-D-O-D-D-U-O-C-C -C-C 
.F-f-f-E-C-t-C-E-£-E-E-E-E-E-E-E-C-E.O-D-D*U-D-U-0-D-C 
.p-F-f -F-F-F-F-F-F-f -F-(-E-£ - E-E -E -L-L-L-E-E -D-O'U-D-O' 

.F-F-F-f-F-f-F-F-F-F.F-F-F-F-F-F-F-E-E-E.C-L.L-E-E-D-D 

»B.li-6-6-6-6-»-6.F-F-f-F-F-F-F.F-F-F-l-F-F-F-E-E-E-E-E' 
-&-a-6-U-6-6-*-6-U-6-6-6-6-f;-6-&-6-6-F-F-F-F«F-F-F*F-E' 
.H-H-H-l*-l*-6-6-6»6-6-6-6-6-6-6-6«6-6-6.6-&-6- 6-F-F-F  «F 

-l1-M.M-M-H-ll-l1-l(-l1-H-H«l1-M-ll-M.M-H.H.H-M»l-l»-&-6-6-6-6 

■  1-ii-*-II-|I*|I-M*M-H»A-H-II*II-*I-iI— H*M— ll.l|*H,ll,H*H*N,H*»*8 
.J.] -1-1 -I-l.|.l-I-l-l-l-l-l-l-I-I-l-l-l-I-I-l-l1-H-H-M 

-1- 1-1- 1- l-l-I-l-l-l-I-l-l-I-l-l- l-l-I-I-l-I.I-I-I-l-I 
■-a-d-d-a-a-a-j-d-a-a-a-a-j-a-a-a-a-a-a-a-j-d-d-a-a-a- l 
-a-a-a-M-K-A -r-r-r-r-m-a.r-k-r-a-r.m.r.a-r-r-r-k-a-r 


1 . 

it  a  e  .  .  . 

MIME*. 

H  H  II  II  f  . 

6  6  H  6  H  o 
6  6  6  6  F  6 
F  6  F  F  C  H 
F  f  f  C  E  H 
F  E  C  C  0  6 
E  L  E  0  C  E 
[EDIICO 
0  0  0  U  C  0 
U  U  l  L  C  C 

0  c  c  c  c  c 

C  C  C  C  B  E 
C  C  C  C  U  c 
CClbM 
C  D  b  U  b  b 
B  B  B  0  b  b 
B  B  b  B  b  B 
b  B  B  B  ■  B 
B  0  B  U  B  n 
8  b  b  b  u  b 

B  B  A  0  b  b 
A  A  A  A  A  B 
A  A  A  A  A  A 
A  A  A  A  A  A 
A  A  A  A  A  A 
A  A  A  A  A  A 
•••AAA 
.A-A  .  .  A  A 
.A-A- A- A  -  A 
a-a-a-a-a  - 
b-A-A-A-A-A 
8-B-U-A-A-A. 
b-b-B-b-b-A. 
-C-C-b-U-b-B- 
-C-C-E-U-b-b- 

-o-c-c-e-c-b- 

l-D-D-D-C-C-C. 
-t -O-U-U-C-C  < 
-E-E -E -D-D-O" 
-F -L-C -L-E  *  O' 
-F-F-F -F-ll-C' 
-6-m-F-F -F-F- 
■H*6 -G-G-b -F 
-H- 11-11-11-6-6' 
1-1-1-I-M- 

-a-j-a-1 -1- I 


E  M 
E  H  II 
0  F  H 
0  E  H 
C  E  6 
CDF 
C  D  E 
C  E  E 
BED 
BCD 
b  C  0 
BBC 
BBC 
BBC 
U  B  B 
ABB 
ABB 
A  A  b 
A  A  b 
A  A  • 
AAA 
•  A  A 
-A-A 
-A-A  A 
A-A. A 
B-A-A 
U-B-A* 
■B-B-b. 

■  C  -B-b* 
■C-C-C* 
•0-O-C* 

■  E  -0-0' 
-l-E-L 
-F-F-E 
i-6-F-F 
l-H-6-6 


-A -A -A -K-A-A-A-K-R-A-K-K-i 


L-L-L-L-L-L-L-L-A-R-K-K-K-a-a-a- 1-1-1 


II  H 
6  K 
F  H 
L  H 
L  6  6 
U  F  6 
DEC 
D  E  6  6 
C  D  F  F 
C  D  E  F 
C  D  E  F 
C  C  0  F 
0  C  0  C 

b  c  c  n 

B  0  C  D 
U  b  C  0 
A  0  C  C 
A  0  B  C 
A  A  B  C 
•  ABB 
-A  A  A  A 
-A-A  A  A 
•b« A-A  A 
-B-b-A-A- 
C-b-B-B- 
■C-C-C-C- 

u-u-o-o- 

-t-E-E-C- 

-F-F-F-F- 

b-6-6*F" 

H-H-H-M. 

1-l-l-H- 

i-a-a-a-i 

a-a-a-a 

-K-j-a 

l-M-K-R 

i-O-b-b 

-R-R-R 

I-S 


0  0 
c  C 
a  c 
b  c 

A  B 
A  A 
'A-A 
•b-B 

•c-e 

•u-c 

•0-D 

•E-E 


I-R-R-R-* 
l-L-K-K-A' 
50  1 -L-L-L-L  < 

1 -m.m-m-l • 

I -H-H-M-R- 

i-h-r-n-r* 

1 -R-K-N-N* 
15  1 -O-O-O-O  ■ 
1-0-0-U-U' 
J-p-p.p-p. 

i-p-p-p-p- 
1-y-B-w-U' 
10  1 -R-O-W-M' 

1-R-fi-R-R' 

1-S-6-N-V 
I-5-SI-S-6 
I-l -6-1- 1 ' 
35  I-f-f-I-1 
I-U-U-U-U' 
1 -U-U-U-u 
I  - v-V-V-V 


A-A.R- 


i-a-a-r-a-r- 


-R-R-R-R-A- 

l-l-l- 

M-H-M- 
M-M-M- 
M-N-K- 
■b-R-M" 
,0-U-O- 
0-0-0- 
F-P-F*' 

o-y-o  < 

O-U-U- 
A-N-K' 
R-R-K- 

s-s-s 

6-S-6' 
1-T-T 
1-T-U' 
U-U-U' 


M-R-1 

R-R-l 

O-U-i 

U-U-l 


■A  —14  —  I 

•R-R-l 

b-Si- 

i-b-i 

T-I- 


m-m-m-m.m-m- 

■R-N-N-N-R-M" 

■M-b-R-R-O-0' 

o-o-o-o-o-a- 

p-p.p-F-p.p. 

.p.p.p-p-p-p 

•w-o-a-o-o-M 

•y-j-o-a-R-A' 
-R-R-R-R-R-R 
■M-R-S-S-6-b 
•S-S-S-S-S-b 
■S-S-I-l-T-l 
■1-T-l-T-l-U 
•U-U-U- U-U-U' 
■u-u-v-v-v-v 
■v-v-v-v-v-y 


•L-L-L-L 

■L-L-L-L* 

■M-H-M-M- 

-M-M-R-R* 

.p-b-b-b* 

■o-o-o-o- 

O-D-f’-P- 

■p.p.p.p. 

•o-y-o-a- 

•O-U-U-R' 
.r-R-R.R- 
•A -S-3-S- 
S-b-S-T 
T-t-t-T 
-T-l -U-U 
•U-U-U-V' 
v-V-V-V 


L-L-L-L-L-L-L-L-t-L-L-L-L-M-M-L-L-L-L-L- 
-M.r.M.M-M-M-M.M-M-M.H-M.M.H-M.M-M. 
-H-N-H.b-M-N-R-R-b-N-b-R-N-b-N-H-b- 
-H-b-R-N-b-O-O-D-O-O.O-O-O-O.O-U-O- 

-o-n-O-D-D-O-P-P-P-F'-P-P-w-B-a-u-F'* 

-p-p.p-p-p-p-p-y-o-y-u-6-R-fl-K-K-K* 

-P-P-U-b-U-W-U-R-A-R-M-M-b-b-i-6-b- 
-0-0-0-R-R-H-M-P-S-S-6-6-T-T-I-T-1- 
-r-r-a-r-S-A-S-6-T -T-t-T-U-U-U-U-U' 
-S-S-S-S-6-l-T-T-T-U-U-U-U-V-V-V-V- 
—  S  —  ft— T  —  1  —  T  —  1 -U-U-U-U- V-V -V-V-M- T - 1 ' 
-T-T-T-U-U-U-U-V-v-V-*-*- ***-2-2-2 
-U-U-U-U- V-V-V-V-y-p-A-T-A- l'l 
-U-V-V,V*W-6-M-M-l»-T -2-2-2 
-V-V-V-m-B-m-A- >-2-2-2 
-«-H-b-«-)l -2-2-1-2 

-•- A-K-2-2-2-2 
-2-2 


L-K-K-l 

M-L-L- 

b-R-M. 

D-O-U- 

P-P-P. 

R-tt-U' 

■  6-S-R' 
T-l-T 
•U-  l-M 
■A-A-R 
•»-T 


R-R-d-a- 
L-L 
R-M-M-L' 
R-R-R-M. 
■P.P-D-0 
'U-U- V-P 
■R-R-l* 
■U-U-U 
-•-V 


30  I-A-i 


-  V- V- V- 
-A-A-*-' 


I-*-T-a-a-»-»-< 


110 


SPECIFIC  IN1ENNAL  lnCHGT 


crtLL« 


THC>4  .OOt)UU£-QE  3ECONUS 


SYMOUL 

NAXjNUK  VALUE 
S1MUOL 

MAXIMUM  VALUE 
5  V  MU JL 

MAAIKUK  value 


V  ■  V  UL  *  U  V  H.VQ£*UV  7.J0EA09  9,HE«U3  l,?tl««10  J.IHUi)  i.AM£*IO  l.V7£«10 


i.HUf «1Q 
1 

**.»UE*iO 


2  •  o**£  *  1 0  a.*ME«IU  i.I3£»10  ).3A£«I0  i.EUt«lQ  l.V^LilO  H.0V£«1U  R.32£«I0 

u  V  V  *  f  t 

3.nvL«iO  3.aJL«io  ».3SIE*IU  s.TvEMO  4 ,  o  ue  ♦  i  o  A.aHE«io 


i  A  A  A  A  A 


:  c  E  l  o  u  p 

I  P  U  P  P  O  N 
I  N  M  N  0  N  N 
I  N  M  N  N  N  N 
I  N  M  N  N  N  N 
I  N  M  N  N  N  N 
I  H  M  N  N  H  N 


l  H  H  * 
I  N  N  » 
I  N  0  A 
I  UO  i 
It  IDt  I 
I  0  0  ( 
I  N  o  ( 
Inn; 
1  M  N  t 
TS  l  M  M  I 


\  N  H 

[  N  N 
»  0  0 
,  0  0 
I  »  0 

)  N  0 
I  N  N 
i  N  N 


N  N  N  N  U  0 
N  N  H  N  0  a 
n  n  a  o  o  a 
a  n  a  u  u  a 
n  n  o  o  u  a 
n  n  a  o  u  a 
n  n  a  o  a  n 

N  N  U  U  N  N 
N  N  a  N  N  N 
N  N  N  N  N  M 

P  M  N  N  M  M 


a  p  u 
o  a  n 

d  N  N 
N  N  N 
N  N  U 

n  n  a 
n  n  a 
N  0  u 
BOO 
a  a  j 
a  o  o 
u  P  B 
b  p  a 
o  b  a 
n  a  a 

N  N  N 
N  N  N 


A  A 

A  A  A  A  A 
A  A  A  A  A 

A  A  A  C  P 

C  H  K  w  a 
g  P  p  p  p 
P  U  P  P  P 

o  a  u  p  p 
a  o  a  p  p 
a  u  o  p  p 
a  o  a  p  p 

o  o  P  P  P 

o  u  p  p  a 
a  a  p  p  p 
a  p  p  p  p 
a  p  p  p  p 
o  p  p  P  P 
a  o  a  u  o 
o  u  u  o  u 
o  a  o  o  a 
o  a  u  o  o 

N  N  N  N  N 
N  N  N  N  N 
N  M  N  N  H 
M  M  N  N  N 


AAA 
AAA 
A  A  0 
CMS 
M  ft  ft 

p  a  m 
u  u  a 
p  a  a 
p  o  a 
u  a  u 
a  a  m 
j  a  a 
n  a  a 

O  U  0 

o  a  a 

p  p  p 
p  p  p 
PUP 

a  a  a 
u  a  o 
u  a  a 
a  a  a 
n  a  a 


AAA 

AAA 

a  a  a 
AST 

p  s  s 

PS# 

ass 
M  S  S 
K  s  S 


W  a  p 
a  o  o 
a  a  o 
p  p  p 

p  p  p 
p  p  p 

o  P  P 

o  a  o 
a  o  a 
a  a  a 
a  a  a 
n  o  o 
n  n  a 

N  N  M 


.  N  M  L  I 


A  A 
A  A  A 
A  A  D 

a  e  o 

b  o  V 
SUV 
T  U  U 
T  U  U 
T  T  u 
s  T  T 
S  T  1 
s  s  s 
s  s  s 
s  s  s 


OOP 
u  a  o 
guv 
p  p  g 
p  p  p 
p  p  p 
p  p  p 

o  p  p 

OOP 

a  a  a 
a  a  u 
a  o  o 
a  a  a 

N  N  0 
N  N  N 
N  N  N 
N  N  N 
N  N  N 


•  *  A  . 

-  A  t  U 

A  |)  V  V 

B  U  A  T 

a  v  x  x 

V  1  V  V 

V  U  V  V 

V  U  V  V 

U  U  v  V 

u  u  u  u 

T  U  u  u 

T  U  U  T 

T  T  T  I 

T  T  T  T 

S  I  S  T 

S  S  S  S 

s  s  s  s 


H  K  A  M 
y  M  M  ft 
JUNK 

a  a  a  u 
p  u  a  a 
p  p  a  o 

p  p  p  o 
p  p  p  p 
p  p  p  p 
p  p  p  p 

o  a  p  p 
a  u  a  a 
o  a  a  o 
a  a  a  a 
n  o  o  u 

N  N  N  N 
N  N  N  N 
N  N  N  N 


V  u  •  . 

V  U  U  i 
A  l  A  \ 
A  7  2  1 

V  T  Y  I 


I  L  L 
I  L  L 
I  L  L 
till  L 


I  K  K 
SO  I  K  K 

I  K  K 


so  i  a  a 
i  a  a 
i  j  j 


LMLLMLLLI 

LMLUMLLLl 

kLULULI 

LLLLLLLLl 

LlLKLLLLI 

ElEKKKKKI 

EkEKHKKKI 

lkkkkkkki 


K  K  K 
K  K  K 
K  K  J 
J  j  J 
a  J  J 
J  o  J 


L  l  L  l  m 

.  L  L  L  L  L 

L  L  L  L  L 

.  K  L  L  L  L 

;  K  K  K  K  E 


K  *  K  J  J  J  J 

J  J  a  J  J  J  a 
a  a  a  j  a  a  a 
a  a  a  a  a  a  a 
a  a  a  a  j  a  a 
j  a  a  a  a  a  a 
a  a  a 


i  l  l 
I  I  I 
i  l  I 
l  l  I 


I  I  I 
I  l  l 
I  l  l 


I  I  l 
a  a*  a 
a  a  a 
a  a  a 
a  a  a 
i  a  a  a 


a  a  a 
a  a  a 
a  a  a 
a  a  a 
a  a  a 
a  a  a 


a  a  K  a  a 
a  a  a  a  a 
j  j  j  .i  ^ 
a  a  a  a  a 
a  a  a  a  a 
j  a  a  a  a 
i  a  a  j  I 
i  i  I  i  i 
I  I  1  l  I 
i  1  i  i  1 
i  I  I  l  I 
i  I  I  i  I 
I  I  I  i  I 
i  I  I  i  I 
i  I  1  i  i 
i  I  i  1  i 
i  I  i  i  i 
1  I  I  i  i 
I  I  i  i  i 
i  i  i  i  a 
1  l  I  a  a 
I  a  a  a  a 
a  a  a  a  a 
a  a  a  a  a 
a  a  a  a  a 
a  a  a  a  a 
j  j  j  j 
a  a  a  a  a 
a  K  k  K  k 
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AABBBBBuCBBCCCCCCDEEE 
BAAAUliUBULABCCCOODO 
BAIINBB  BBCCCCC 
A  B  B  A  B  B  B 


A  A  B 

A  A  B 

ABB 
A  8  8 

BBS 

b  El  B 

W  ft  c 

c  c  c 

c  c  c 

C  C  0 

D  D  0 

D  0  U 

0  E  E 

0  E  E 

E  E  F 

F  F  6 


-A-B- 

-A-B- 
-A-C-B- 
-C-C-B" 
-O-C-B- 
■E-O-C-8* 
■0 -O-C-B 
►O-U-C-B 

-E-O-C-A 

-D-O-C  - 
■O-C-B  A 
-O-C-A  a 
■0-B  -  ti 
■C-B  A  C 
■B-A  B  C 
■B  A  B  C 
•A  A  C  0 
-  U  C  0 
a  a  c  e 

A  C  0  E 
8  C  0  E 
B  C  C  F 
C  0  E  F 
C  0  E  F 
C  D  E  F 
0  t  F  C 
0  L  F  6 
0  C  F  6 
E  E  F  6 
E  F  F  6 
£  F  G  6 
E  F  6  & 
E  F  0  H 
C  F  '6  H 
F  F  6  H 
F  F  6  H 
F  I  &  M 
F  F  6  H 
F  F  6  6 
F  F  0  6 
F  F  &  6 
I  F  6  G 
E  F  6  6 
E  F  F  C 
C  F  F  C 
C  F  F  G 
C  F  F  G 
E  E  F  G 
E  E  F  6 
t  t  F  F 
E  E  F  F 
0  E  F  F 
0  E  F  F 
0  E  E  F 
0  E  E  F 
D  E  E  F 
0  0  E  F 
0  0  E  E 
0  0  E  C 
0  0  E  E 
C  0  E  E 
C  0  0  C 
CODE 

coot 

C  0  0  E 
C  C  0  D 
C  C  0  D 
C  C  0  0 
C  C  0  D 
B  C  0  D 
B  C  C  0 

•  CCD 
B  C  C  0 
B  C  C  0 

•  CCD 
B  B  C  C 
B  B  C  C 
B  B  C  C 
G  B  C  C 
B  B  C  C 
B  B  C  C 
B  B  C  C 
B  B  C  C 
B  C  C  C 
8  C  CD 
B  C  C  0 
e  C  C  D 
C  c  0  D 
C  D  0  D 
D  D  D  E 
D  D  C  C 
D  C  E  E 
E  C  E  E 
E  F  F  & 
E  F  C  6 
G  G  H 


A  A  C  C  . 

B  -  B  C  0  I 

■6  A  B  C  E  I 

A  A  B  C  E  . 

A  A  B  C  F 

A  A  B  0  G  . 

A  B  C  C  II  I 

-  0  0  F  K  I 

A  C  D  F  I  I 

BCCGI  I 
e  0  E  G  I  I 

C  0  F  H  J  I 

C  E  F  H  J  I 

D  E  G  H  J  I 

0  E  6  I  J  I 

D  F  C  I  J  I 

E  F  M  I  K  I 

EiHlKI 
F  G  H  J  K  I 

F  6  H  J  K  I 

F  6  I  J  K  I 


H  1  J  K  L  I 

H  1  J  K  K  I 

M  I  J  K  K  I 

H  1  J  K  K  I 

H  1  J  K  K  I 

H  1  J  J  K  { 

H  1  J  J  K  I 

Hi  J  J  K  I 

HI  J  J  K  I 

H  I  I  J  K  | 

K  I  I  J  K  I 

Hi  I  J  K  I 

HI  I  J  K  I 

H  M  I  J  K  I 

H  li  I  J  J  I 

H  H  1  J  J  I 

H  II  I  J  J  I 

6  K  I  I  J  I 

G  K  I  I  J  l 

G  li  I  I  j  I 

G  H  H  I  J  . 

€  II  H  I  J  , 

G  II  K  I  J  , 

G  G  H  I  I  , 

G  G  II  1  I  , 

G  6  H  I  J  . 

F  G  H  H  I  , 

F  G  II  II  I 

F  6  G  II  I 

F  G  G  H  I 

F  6  0  li  K 

F  F  G  II  H 

F  F  6  H  II 

F  F  6  6  h 

E  F  G  G  H 

CFCGHI 
E  F  F  G  II  I 

E  F  F  G  6  I 

C  F  F  G  G  I 

E  E  F  C  G  I 

E  E  F  F  6  I 

E  E  F  F  6  I 

0  E  F  F  G 

D  E  E  F  G  I 

0  E  E  F  G  ' 

0  E  E  F  F 

0  E  C  F  F 

D  0  C  F  F  I 

D  0  E  E  F 

0  0  E  E  F 

0  0  C  E  F  I 

0  D  E  C  F  I 

D  0  C  E  F  I 

0  D  C  E  F  i 

0  0  C  E  F  I 

0  0  E  E  F  i 

0  0  E  F  Fi 

0  E  E  F  F  i 

0  E  E  F  G  I 

C  E  F  F  6  i 

E  E  F  F  G 

E  F  F  G  FI 

E  F  &  6  & 

F  F  G  M  H 

F  F  li  H 

6  0  H 


>  P  U 

:  o  R 

INST 
i  N  R  t 

1  H  0  U  U 

l  H  P  I  V 

.  U  P  6  v  l 

.  *  P  *  U  i 

.  H  P  R  T  i 

.  U  P  0  S  i 

.  h  0  0  S  I 

.MOOS 
.NOOK 
.  M  0  0  K 

.  M  o  r  fi 

,  H  0  P  H 

.  N  0  P  R 

■  H  0  P  0 

I  H  0  P  U 

I  H  0  P  U 

.  H  D  P  U 

.  H  0  P  0 

.  M  M  O  P 

.  M  N  U  P 

.  M  H  0  P 

.  M  H  0  P 

.  M  N  0  P 

.  M  N  0  P 

.  ¥  H  0  P 

.  M  H  0  0 

.  M  H  H  0 

.  K  M  H  U 

.  L  K  H  0 

.  L  H  H  0 

.  L  *  H  D 

I  L  M  H  D 

<  L  K  N  H 


,  L  H  H  M 

L  L  M  M 
L  L  M  M 
K  L  H  N 
K  L  M  H 
K  L  M  M 
'  K  L  H  M 

'  K  L  L  M 

I  K  L  L  K 

'  K  K  L  M 

'  K  K  L  M 

i  J  K  L  M 

I  J  K  L  L 

J  K  L  L 


I  J  K  K 
I  J  K  K 
I  J  J  K 
I  J  J  K 
I  I  J  K 
I  I  0  K 
H  I  J  J 
H  I  J  J 
H  I  I  0 
H  I  I  J 
H  li  I  J 
H  H  I  J 
H  H  I  I 
G  H  I  I 
6  H  I  I 
G  H  H  I 
G  H  H  I 
G  H  H  I 
G  H  H  I 
G  H  H  I 
0  H  H  I 
G  H  I  I 
6  H  1  I 
H  H  I  J 
H  I  1  J 
H  I  J  J 
H  I  J  I 
H  I  J  I 
H  J  K 


I  V  7  Y 

I  U  M  Y 

S  U  V  T  2 

S  I  V  X  2 

S  T  V  M  2  2 

X  T  u  V  *  2 

H  I  U  V  -  2  3 

H  5  U  V  «  A  i 

K  G  I  U  h  X  ; 

M  G  I  U  V  G  j 

K  S  T  U  V  h  1 

w  R  S  u  V  w  I 

w  H  S  T  u  V  I 

0  R  S  T  U  V  « 

0  R  S  T  U  V  l 

y  R  s  I  U  v  i 

P  0  ft  G  T  U  \ 

P  0  K  S  I  U  I 

P  0  R  S  I  U  \ 

P  0  R  S  T  U  \ 

P  U  K  S  S  T  I 

P  0  0  M  $  1  I 

p  P  0  R  S  I  l 

DP  t  II  $  I  I 

0  P  0  R  $  T  i 

0  P  0  K  S  6  1 

0  P  U  R  R  S  1 

D  P  0  0  R  6  1 

0  P  p  o  H  S  ' 

U  0  P  Q  R  S  ’ 

0  D  P  0  R  S  ! 

H  0  P  0  K  H  < 

H  0  P  o  0  K 

M  0  P  R  «  M  J 

M  0  0  P  Q  R  I 

M  0  0  P  0  H  ! 

M  U  0  P  Q  R  I 

K  W  0  P  0  b  I 

H  H  0  P  P  0  I 

M  H  0  P  P  b  I 

K  M  0  0  P  U  I 

M  M  M  0  P  0  I 

M  M  M  0  P  0  ( 

M  M  H  0  P  P  ( 

L  M  H  0  0  P  I 

L  M  M  H  D  P  I 

L  M  N  H  0  P  I 

L  H  M  M  D  P  I 

L  L  M  k  0  P  I 

L  L  M  H  0  0  I 

K  L  M  M  M  0 

K  L  H  m  M  0  I 

K  L  L  M  M  0  I 

K  L  L  M  h  0  I 


J  K  l 
J  J  I 
J  J  I 
J  J  I 
J  0  I 
J  J  l 
J  J  I 


K  L  L  I 
K  K  L  I 

L  L  K  i 
L  K  K  I 
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AXIAL  VC  LUC  ITT 


CYClC*  I  tt7 . 0  TIMf*5.21lYIk.-t>&  SCCUND6 


fYMdUL 

MAXIMUM  VALUC 
SYMBOL 

MAXIMUM  VALUC 
fiYNHUL 

MAXIMUM  VALUE 


l.»hC»U3  I.5ALY0H  l,m«OH 


•  ctc*  dh  g.21e*qh  t.yi  l«oh  h.I5E»uh  i.qic«05  i.^sefos 


I.iUL»05  I.56E*t»5  I.VILMJS  I.v7t*04  a.0JC»05  2.I«L*US  S.JHC«05  2.H»t»U3  2.6iE«05  2.#0C«05 

6  1  II  V  w  X  Y  i 

Z.16E*QS  A.I2C*V3  1.27L405  l.nlC*05  1.5»£*05  i.TIE^OS  5.1QE»0S  h.u5l*US 


100  I  M  M  M  M  N  M  I 

I  M  M  M  M  M  L  I 

1  L  L  L  L  L  L  I 

I  R  R  R  R  K  R  I 

I  K  K  K  A  K  R  I 

15  I  d  J  J  J  J  J 

I  1  I  I  I  I  I 

I  H  M  H  K  K  H  I 

I  H  H  H  H  «  * 

I  U  G  A  6  6  5 

»0  I  f  F  F  f  F  F 

If  FFf  FF 
I  C  C  C  L  L  C 
I  E  E  L  C  t  U 
I  L)  U  U  L>  U  U 

US  I  U  o  U  u  c  c 

ICtCckC 

I  b  B  b  d  b  b 
IMbUIB 

I  A  A  A  A  A  A 

•0  I  A  A  A  A  A  A 

I -A -A- A- A-A- A < 

I -A-A-a-A-A-A- 

| -O-B-d-b-b-U- 
J-B-b-O-L-C-C. 
75  I-C-C-C-l-L-t- 
l-O-O-U-U-U-U- 
I-O-D-O-O-O-l- 
l-E-E-E-l-L-L- 
I-E-E-C-L-L-E- 
TO  I -F -F-F-F  -F-F. 
I- f-f -f-f-f-f- 

j-F-F-F -F-F-F. 
j-F-F-F -F-F-F. 

1  -6  -6-6-6 -6-l>« 
65  I  -6-6-6-6-6-G. 

j-6-e-ti-t.-i.-6. 

I-6-G-6-G-B-6. 

I -H-H-6-6-M-6* 
J-H-H-H-II-M-H. 
60  I  — H— M— H— H— H-H< 

I  — H— M— H— II— H-ll* 
I-I-I-I-l-I-l. 
I-I-I-I-I-l-I 
I-l-I-I-I-l-I 
55  I-J-d-J-J-d-d. 

I -d-d-J-d-d-d' 
J-J-d-J-J-J-d 
I -k-k-m-k-m-k. 


J  J  J 
I  I  I 
H  H  II 
H  H  6 
G  6  6 
6  F  F 
FFF 
L  l  t 
L  C  C 
U  U  Li 

0  0  c 
C  C  L 
CCU 


c-c-c- 

c-c-c 

o-o-u. 

>0  -U-U- 

l-c-c. 

■E-E-L . 

F-F-F. 

F-F-F 

■F-F-F. 

■F-F-F. I 

■F-F-F-l 

■G-6-6-1 

■  U-H-b.l 
■6-to-b-i 
■6-G-U- 
•H-H-H-l 
.H-ll -I  I -I 
l-IF-ll-tl-l 
-I-I-l- 
-l-I-I- 
-  I  -  I  - 1  - 
-d-d-d- 
-d-d»d- 
l-d-d-d- 


L  K  *  IU  J 

R  R  R  J  J  J 
A  d  J  J  J  I 
J  d  d  I  1  H 
J  1  I  H  H  6 
I  K  K  H  6  6 
H  R  6  G  G  F 
G  6  G  F  F  F 
6  F  F  F  F  C 
F  F  F  E  C  C 
feme 

C  c  C  C  D  D 
C  0  0  0  U  0 
D  o  0  U  C  C 
0  c  c  c  c  c 

C  C  C  C  C  B 
C  d  B  b  8  e 
b  u  b  v  n  a 


H  H 
I  I  H  H 

I  I  H  H  6 

J  I  !  H  H  6 

J  I  HHH( 

I  I  h  e  6  f 

I  H  6  6  F  F 
H  6  %  F  F  E 
6  6  F  F  E  E 

6  F  F  C  C  C 

F  F  C  E  £  E 

F  E  C  E  C  0 

E  C  E  C  0  0 

E  C  £  0  0  0 

c  o  o  n  o  o 

0  D  0  D  D  C 

0  D  D  c  C  C 

o  c  c  c  c  c 

C  C  C  C  B  B 

c  c  B  B  B  b 

U  0  D  B  6  B 

B  B  R  B  A  A 


C  0 
CEO 
F  F  E  E 

G  F  F  £  E 

6  6  F  F  E  E 

6  6  F  F  E  D 

»  F  F  E  E  0 

6  6  F  E  £  D 

6  F  E  E  0  U 

F  F  E  E  0  0 

F  E  E  E  D  0 

E  E  E  0  0  0 

£  E  C  U  0  U 

£  £  E  0  D  0 

E  E  0  D  0  0 

C  0  0  D  o  0 

D  0  0  0  0  C 

0  D  0  0  C  C 

0  0  0  C  C  C 

DCCCCt 
C  c  C  C  B  B 
C  C  C  H  S  B 
C  C  B  B  B  B 
II«B(  A 
6  fi  b  A  A  A 
B  B  A  A  A  A 
A  A  A  A  A  - 
AAA  -A-A 


A  C  6  6 

B  0  6  6 

C  0  F  6 

C  F  F  6 

C  F  6  F 

0  E  f  E 

E  C  F  E 

E  E  F  E 

DDCDEECD 
0  0  0  E  C  C  0 

0  U  D  E  £  D  D 

0  0  D  c  C  0  D 

D  0  E  E  E  D  D 

0  D  U  E  0  0  C 

0  0  D  D  D  D  C 

D  O  0  D  D  C  C 

0  0  0  0  0  c  c 

0  0  D  0  0  C  C 

0  D  0  0  C  C  C 

0  D  0  C  C  C  c 

0  D  0  C.  C  C  B 

0  0  C  C  C  C  B 


i  i  u 
!  B  b 

I  B  U 
I  S  fi 


1  I  E  . 

H  6  6  E 

F  F  F  6 

f  E  L  H 

E  0  D  K  I 

0  D  C  F  I 

D  D  C  0  I  d 

D  0  C  D  6  I 

0  C  C  C  E  I 

D  C  C  C  E  H 

0  C  C  C  0  F 

C  C  C  C  D  E 

C  C  C  C  D  E 

C  C  C  C  C  o 

c  c  c  c  c  c 

C  C  C  C  C  D 

c  c  c  c  c  c 

C  C  «  B  C  C 

C  B  B  b  C  C 

B  B  b  •  •  C 

B  B  b  B  B  C 

b  B  t>  H  b  b 

B  B  a  b  B  • 

J  B  W  b  4  B 

b  B  d  Q  li  8 

d  b  B  B  B  B 

A  A  A  B  B  b 

A  A  A  A  B  B 


A  A  A  A  A  i 


•  •  A  A  A  A  A  A  A 


--A-A- A-A-A 
- A-A-A- A-A-A-i 


d  d 
H  d 
F  d  « 

F  H  J 
E  6  J  d 
E  F  R  d 
D  f  6  d 
0  E  F  IF 
D  t  F  6 
C  0  E  F 
C  D  E  F 
C  U  D  E 
C  C  D  E 
CEDE 
•  i  D  0 
B  C  C  D 
HCD 
b  u  c  c 
b  B  C  C 
b  B  B  C 
S  b  •  C 

d  b  a  c 

A  w  B  B 
A  B  •  b 
A  b  B  B 
A  A  b  B 
A  A  B  B 
A  A  B  B 


A  A  A  A  A  - 
-  —• A-A-A— A< 
.A-A-X-A-A-A' 
.A-b-n-b-u-M 
■B-D-H-d-d-B' 
C-L-C-L-C-C 
C-C-C-C-C-C 

■rt-o-n-u-o-ti- 
•o-u-n-u-u-o. 
n-u-n-E-o-o- 
C-l-C-E-C-E 
■t  -l-c  -l  -c -E  ' 
C-L-F-L-L-C. 
■F-F-F-L-E  -£< 
F-F-F-F-f-E 
F-F-F -F-F-F 
■F  -I  -F-F-F-F 
F-F-F-F-F-F 
6-G-G-G-F-F 
G-6-6-6-G-G 
■C-6-G-6-G-G 
>11-6-6-6-6-6 
■iuh-h-h-m-i 

•  R-H-H-H-IF-I 
■  I — I  —  I  -rt-H-i 
•I-I-I- I-I-l 
-I-l-I-I -I-I 
>d-d-d-d-J-d 
•J-d-J-d-d-J' 

•  d-d-.l-d-d-J' 


,-A-A-A-i 

-a-b-b-1 

B-M-b. 

B-8-17. 

c-c-c- 

c-c-c. 

C-C-C. 

D-n-c. 

F)-D-0 . 

D-n-D. 

o-d-d- 

■O-D-O. 

E-n*D. 

E-E-E- 

E-r-E 

E-F-E 

F-F-C. 

f-f-f. 

■F-F-F. 
G-S-f . 
•6-G-6 


i-G-t 


-h-k-h-h-H' 


h-h-b 

-H-IF-H-R-R-H-H. 

ii-h-if. 
I-I-I. 
1-7  -  I 
d-d- l 


d-d* 

-d-d- 

d-d- 


|-d- 

d-.l-J- 


a-b-b- 
b-B-' ■ 
B-B-B. 
H-B-d. 
C-C-C- 
•C-C-C. 

c-c-c. 

C-C-C- 

C-C-C. 

o-o-c. 

D-D-O- 

U-D-O. 

O-U.D. 

E-O-U- 

E-E-E. 

E-t-C. 

E-E-E- 

•  F-F-F- 
F-F-F- 
•F-F-F- 

6-6-6- 

6-6-6. 

6-6-6 

■H-K-H< 

■M-M-M. 

•  I-I-K. 
I-I-I 
I-I-I 
d-d-d- 
J-d-d 


B-B-B-B 
B-b-B-B 
■B-d-B-d 
.U-b-B-D 
C-B-B-B 
C-C-B-U- 
■  C.C-C-B 

•C-C-C-l 
.c -c-c-c 

C-C-C-C 

o-c-c-c 

u-o-o-c 

U-U-O-D 
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Appendix  C-IVa  Copper  Wedge  Startup  Output 

Here,  the  compression  map  is  shown  at  cycle  zero  to  show  the  initial 
configuration  of  the  copper  wedge  at  impact  time.  At  the  time  p/pq  =  1 

for  all  cells.  A  compression  map  is  also  shown  at  1.9  usee  (cycle  11). 
Here  some  compression  can  be  seen  at  the  reflecting  wall  boundary. 
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COPPER  IMPACT  STUDY  6Q  DEGREE  WEDGE 
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Appendix  C-IVb.  Copper  Wedge  Restart  Output 


At  .295  ysec  (cycle  14)  compressions  can  be  seen  to  be  increasing 
at  the  reflective  wall,  and  the  specific  energy  map,  also  shown  here, 
indicates  specific  internal  energies  in  excess  of  3.73  E+10  ergs/g. 
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Hampton,  VA  23365 

1  Director 

National  Aeronautics  and 
Space  Administration 
Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 


No.  of 

Copies  Organization 

2  Director 

National  Aeronautics  and 
Space  Administration 
George  C.  Marshall  Space 
Flight  Center 
ATTN:  MS- 1/ Lib 

R-AERO-AE/A.  Felix 
Huntsville,  AL  35812 

1  Cal span  Corporation 
P.0.  Box  235 
Buffalo,  NY  14221 

1  Lockheed  Missile  and 
Space  Company 
ATTN:  Mr.  J.  E.  May 
55-80,  Bldg  57 
P.0.  Box  504 
Sunnyvale,  CA  94088 

1  Martin  Marietta  Aerospace 
ATTN:  Mr.  Culotta 
P.0.  Box  5837 
Orlando,  FL  32805 

1  Orlando  Technology  Inc. 
ATTN:  Mr.  D.  Matuska 
P.0.  Box  855 
Shalimar,  FL  32579 

1  Physics  International  Corp. 
ATTN:  Mr.  L.  Behrmann 

2700  Merced  Street 
San  Leandro,  CA  94577 

1  Rockwell  International 
Science  Center 
ATTN:  Dr.  N.  Malmuth 

P.0.  Box  1085 
1000  Oakes,  CA  91360 

1  Science  Applications  Inc. 
ATTN:  Ken  Dent,  Suite  800 
2019  W.  Clinton  Avenue 
Huntsville,  AL  35805 


1  ASD/XRA  (STINFO) 

Wright-Patterson  AFB,  OH  45433 
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DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 

1  Systems,  Science  §  Software 
ATTN:  Dr.  R.  Sedgwick 
P.0.  Box  1620 

La  Jolla,  CA  92037 

2  University  of  California 

Los  Alamos  Scientific  Laboratory 
ATTN:  Dr.  R.  Karpp 
Tech  Lib 
P.0.  Box  1663 
Los  Alamos,  NM  87545 

2  Drexel  Institute  of  Technology 
Wave  Propagation  Research  Center 
ATTN:  Prof.  P.  Chou 

Dr.  J.  Carleone 
32d  5  Chestnut  Streets 
Philadelphia,  PA  19104 

1  Director 

Applied  Physics  Laboratory 
The  Johns  Hopkins  University 
John  Hopkins  Road 
Laurel,  MD  20810 

1  Southwest  Research  Institute 
ATTN:  Mr.  P.  S.  Westine 
P.0.  Drawer  28510 
8500  Culebra  Road 
San  Antonio,  TX  78228 

1  SRI  International 

ATTN:  Dr.  Lynn  Seaman 
Menlo  Park,  CA  94025 

Aberdeen  Proving  Ground 


Dir,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP/H.  Cohen 
Cdr,  USATECOM 

ATTN:  DRSTE-TO-F 
Dir,  Wpns  Sys  Concepts  Team 
Bldg.  E3516,  EA 
ATTN:  DRDAR-ACW 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet  and  return  it  to  Director,  US  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  ATTN:  DRDAR^TSB,  Aberdeen  Proving  Ground, 
Maryland  21005.  Your  comments  will  provide  us  with  information 
for  improving  future  reports. 

1 .  BRL  Report  Number^ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name :  _ _ 

Telephone  Number: 


Organization  Address: 


